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Abstract 

In the past years there has been a revival of hadron spectroscopy. Many interesting new hadron states were discovered experimen- 
tally, some of which do not fit easily into the quark model. This situation motivated a vigorous theoretical activity. This is a rapidly 
evolving field with enormous amount of new experimental information. In the present report we include and discuss data which 
were released very recently. The present review is the first one written from the perspective of QCD sum rules (QCDSR), where 
we present the main steps of concrete calculations and compare the results with other approaches and with experimental data. 
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We are approaching the end of a decade which will be 
remembered as the "decade of the revival of hadron spec- 
troscopy". During these years several e + e~ colliders started to 
operate and produce a large body of experimental information. 
At the same time new data came from the existing p — p collid- 
ers and also from the e — p accelerators. In Table[T]we give a list 
of the new charmonium states observed in these accelerators. 

Table 1 : Charmonium states observed in the last years. 



state 


production mode 


decay mode 


ref. 


X(3872) 


B -> KX(3S72) 




UJ 


X(3915) 


yy -> X(3915) 


J/lf/d) 


m 


Z(3930) 


yy -> Z(3930) 


DD 


□j 


7(3930) 


B -> #7(3930) 


JI4ioj 


[4J 


X(3940) 


e + e- -> J/if/X(3940) 


DD* 


[5J 


7(4008) 


e + e~ -> 7(4008) 


Jlif/nn 


[6J 


Z+(4050) 


B° -> /TZ+(4050) 


Xc\n + 


□ 


7(4140) 


B -> #7(4140) 


Jin 




X(4160) 


e + e- -> //^X(4160) 


D*D* 


[9J 


Z+(4250) 


B° K~Z\ (4250) 




W 


7(4260) 


e + e- -> r/SR 7(4260) 




[10] 


X(4350) 


yy -> X(4350) 


Jin 


UJJ 


7(4360) 


e + e~ -> r/sfi7(4260) 


xj/'nn 


U2J 


Z + (4430) 


B° #-Z + (4430) 


ijf'n + 


UJ 


X(4630) 


e + e- -> y/s«^(4630) 


A + A- 


U4J 


7(4660) 


e + e- ^7/5*^(4660) 


i]/'nn 


[12] 



In what follows we will review and comment all this infor- 
mation ffl llJ^I^HJiJl JTSJIilJi 

MElEilMliflEiB. 
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The study of spectroscopy and the decay properties of the 
heavy flavor mesonic states provides us with useful informa- 
tion about the dynamics of quarks and gluons at the hadronic 
scale. The remarkable progress on the experimental side, with 
various high energy machines has opened up new challenges in 
the theoretical understanding of heavy flavor hadrons. 

1.1. New experiments 

The B-factories, the PEPII at SLAC in the U.S.A., and the 
KEKB at KEK in Japan, were constructed to test the Standard 
Model mechanism for CP violation. However, their most in- 
teresting achievement was to contribute to the field of hadron 
spectroscopy, in particular in the area of charmonium spec- 
troscopy. They are e + e~ colliders operating at a CM energy 
near 10,580 MeV. The BB pairs produced are measured by the 
BaBar (SLAC) and Belle (KEK) collaborations. Charmonium 
states are copiously produced at the B-factories in a variety of 
processes. At the quark level, the b quark decays weakly to a 
c quark accompanied by the emission of a virtual W~ boson. 
Approximately half of the time, the W~ boson materializes as a 
sc pair. Therefore, half of the B meson decays result in a final 
state that contains a cc pair. When these cc pairs are produced 
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close to each other in phase space, they can coalesce to form a 
cc charmonium meson. 




ISR 



Figure 1: Illustration for the initial state radiation (ISR) process. 

The simplest charmonium producing B meson decay is: B — > 
K(cc). Another interesting form to produce charmonium in B- 
factories is directly from the e + e~ collision, when the initial 
state e + or e~ occasionally radiates a high energy y-ray, and the 
e + e~ subsequently annihilate at a corresponding reduced CM 
energy, as ilustrated in Fig.Q] 

When the energy of the radiated y-ray (yisn) is between 4000 
and 5000 MeV, the e + e~ annhilation occurs at CM energies that 
correspond to the range of mass of the charmonium mesons. 
Thus, the initial state radiation (ISR) process can directly pro- 
duce charmonium states with J PC — 1 . 

1.2. New states 

Many states observed by BaBar and Belle collaborations, like 
the X(3872), F(3930), Z(3930), X(3940), F(4008), Z+(4050), 
7(4140), X(4160), Z+(4250), F(4260), F(4360), Z + (4430) and 
F(4660), remain controversial. A common feature of these 
states is that they are seen to decay to final states that contain 
charmed and anticharmed quarks. Since their masses and decay 
modes are not in agreement with the predictions from poten- 
tial models, they are considered as candidates for exotic states. 
By exotic we mean a more complex structure than the sim- 
ple quark-antiquark state, like hybrid, molecular or tetraquark 
states. The idea of unconventional quark structures is quite old 
and despite decades of progress, no exotic meson has been con- 
clusively identified. In particular, those with qq quantum num- 
bers should mix with ordinary mesons and are thus hard to un- 
derstand. Therefore, the observation of these new states is a 
challenge to our understanding of QCD. 

In 2003 a great deal of attention was given to pentaquark 
states, which may be exotic objects l42l 14311 , After many nega- 
tive results the study of these states was abandoned. However, 
in that same year the discovery by BELLE of the unusual char- 
monium state named X prompted a lively and deep discussion 
about the existence of exotic states in the charmonium sector 
(i.e. non pure c-c states), exactly where, up to that moment, 
the calculations based on potential models had worked so well. 

Establishing the existence of these states means already a re- 
markable progress in hadron physics. Moreover it poses the 
new question: what is the structure of these new states? The 
debate on exotic hadron structure was strongly revived during 



the short "pentaquark era". Unfortunately, because of the un- 
certainties on the experimental side, it was very soon aborted. 
However, some interesting ideas were passed along to the sub- 
sequent discussion about the nature of the X, Y and Z states, 
which is still in progress. 

1.3. New quark configurations 

Concerning the structure we can say that there are still at- 
tempts to interprete the new states as c - c. One can pursue this 
approach introducing corrections in the potential, such as quark 
pair creation. This "screened potential" changes the previous 
results, obtained with the unscreened potential and allows to 
understand some of the new data in the c-c approach 14411 . De- 
parting from the c-c assignment, the next option is a system 
composed by four quarks, which can be uncorrelated, form- 
ing a kind of bag, or can be g rouped in diquarks which then 
form a bound system J45, 46, 47]. These configurations are 
called tetraquarks. Alternatively, these four quarks can form 
two mesons which then interact and form a bound state. If the 
mesons contain only one charm quark or antiquark, this con- 
figuration is referred to as a molecule 1 48 , 49, 5^1 . If one of 



the mesons is a charmonium, then the configuration is called 



hadro-charmonium (51 
hybrid charmonium 152 



. Another possible configuration is a 
5311 . In this case, apart from the c-c 
pair, the state contains excitations of the gluon field. In some 
implementations of the hybrid, the excited gluon field is repre- 
sented by a "string" or flux tube, which can oscillate in normal 
modes |50|. 

The configurations mentioned in the previous paragraph are 
quite different and are governed by different dynamics. In 
quarkonia states the quarks have a short range interaction domi- 
nated by one gluon exchange and a long range non-perturbative 
confining interaction, which is often parametrized by a linear 
attractive potential. In tetraquarks besides these two types of in- 
terations, we may have a diquark-antidiquark interaction, which 
is not very well known. In molecules and hadro-charmonium 
the interaction occurs through meson exchange. Finally, in 
some models inspired by lattice QCD results, there is a flux 
tube formation between color charges and also string junctions. 
With these building blocks one can construct very complicated 
"stringy" combinations of quarks and antiquarks and their in- 
teractions follow the rules of string fusion and/or recombination 
lf54l l55ll . In principle, the knowledge of the interaction should 
be enough to determine the spatial configuration of the system. 
In practice, this is only feasible in simple cases, such as the 
charmonium in the non-relativistic approach, where having the 
potential one can solve the Schrodinger equation and determine 
the wave function. In other approaches the spatial configuration 
must be guessed and it may play a crucial role in the production 
and decay of these states. 

1.4. New mechanisms of decay and production 

The theoretical study of production and decays of the new 
quarkonium states is still in the beginning. Some of the states 
are quite narrow and this is difficult to understand in some of 
the theoretical treatments, as for example, in the molecular ap- 
proach. 
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The narrowness of a state might be, among other possibili- 
ties, a consequence of its spatial configuration. Assuming, for 
example, that the X(3872) is a set of cc uu quarks confined in 
a bag but with sufficient (or nearly sufficient) energy to decay 
into the two mesons D and D , why is it so difficult for them to 
do so? What mechanism hinders this seemingly simple coales- 
cence process, which is observed in other hadronic reactions? 
In a not so distant past the same questions was asked in the 
case of the + pentaquark. Now we know that the existence 
of this particle is, to say the least, not very likely. Neverthe- 
less, as mentioned in a previous subsection, some of the ideas 
advanced in the pentaquark years might be retaken now in a 
different context to answer the question raised above. In fact, 
the pentaquark decay & + (uudds) — > n{udd) + K + (us) was en- 
ergetically "superallowed", since there was enough phase space 
available and no process of quark pair creation nor annihilation. 
Why was this quark rearrangement and hadronization difficult? 
It has been conjectured 115 611 that the + was in a diquark (ud) 
- diquark (ud) - antidiquark (s) configuration, such that, due to 
the not very well understood repulsive and/or attractive diquark 
interactions, the two diquarks were far apart from each other 
with the antiquark standing in the middle. In this configuration 
it was difficult for one diquark to capture the missing d quark (to 
form a neutron), which was very far. These words were imple- 
mented in a quantitative model and this configuration was bap- 
tized "peanut", since the three bodies were nearly aligned. An- 
other pentaquark spatial configuration with small decay width 
was the equilateral tetrahedron with fourquarks located in the 
corners and the antiquark in the center 157] . An other interest- 
ing conjecture based on lattice arguments was that, due to the 
dynamics of string rearrangement, the pentaquark constructed 
with three string junctions had, during the decay process, to 
pass through a very energetic configuration. Since its initial 
energy would be much less than that, this passage would have 
to proceed through tunneling and would therefore be strongly 
suppressed |54j|. 

These ideas are relevant for the new charmonium physics be- 
cause if they are multiquark states, as it seems to be the case, 
their spatial configuration will play a more important role both 
in production and decay. 

The production of some of these states (such as, for exam- 
ple, the X) has been observed both in e + e~ and p p colliders, 
the latter being much more energetic. The X produced in elec- 
tron - positron collisions comes from B decays, whereas in p p 
reactions it must come from a hard gluon splitting into a cc 
pair, which then undergoes some c omp licated fragmentation 
process. In theoretical simulations 115 811 . it was shown that it 
is very difficult to produce X if it is composed by two bound 
mesons. 

Some reviews about these states can be found in the literature 

lissmmsiiiiisiiiiii]. There 31-6 at least two rea_ 

sons for writting a report on the subject. The first one is because 
this is a rapidly evolving field with enormous amount of new 
experimental information coming from the analysis of BELLE 
and BABAR accumulated data and also from CLEO and BES 
which are still running and producing new data. New data are 
also expected to come from the LHCb, which will start to op- 



erate soon and will be generating data with very high statistics 
for the next several years. In the present report we include and 
discuss data which were not yet available to the previous re- 
viewers. 

The other reason which motivates us to write this report is 
that, on the theory side each review is biased and naturally em- 
phasizes the approach followed by the authors. Thus, in Refs. 
j60ll and l66ll the authors present the available data and then 
discuss their theoretical interpretations making a sketch of the 
existing theories. In IHoll attention is given to the conventional 
quark antiquark potential model and to models of the interac- 
tion between mesons which form a molecular state. In l66ll 
a very nice overview of different theoretical approaches such 
as potential models, QCD sum rules and lattice QCD was pre- 
sented. The author works out some pedagogical examples, us- 
ing general considerations and simplified assumptions. The 
present review is the first one written from the perspective of 
QCD sum rules, where we present the main steps of concrete 
calculations and compare the results with other approaches and 
with experimental data. 

In the next sections we discuss the experimental data and the 
possible interpretations for the recently observed X, Y and Z 
mesons. 

2. QCD sum rules 

2.1. Correlation functions 

QCD sum rules are discussed in many reviews 
7211 emphasizing various aspects of the method. The basic idea 
of the QCD sum rule formalism is to approach the bound state 
problem in QCD from the asymptotic freedom side, i.e., to 
start at short distances, where the quark-gluon dynamics is es- 
sentially perturbative, and move to larger distances where the 
hadronic states are formed, including non-perturbative effects 
"step by step", and using some approximate procedure to ex- 
tract information on hadronic properties. 

The QCD sum rule calculations are based on the correlator of 
two hadronic currents. A generic two-point correlation function 
is given by 



H(q) = if d 4 xe k '-\0\T[j(x)f(0)]\0) , 



(1) 



where j(x) is a current with the quantum numbers of the hadron 
we want to study. 

The fundamental hypothesis of the sum rules approach is the 
assumption that there is an interval in q over which the corre- 
lation function may be equivalently described at both the quark 
and the hadron levels. The former is called QCD or OPE side 
and the latter is called phenomenological side. By matching 
the two sides of the sum rule we obtain information on hadron 
properties. 

2.2. The OPE side 

At the quark level the complex structure of the QCD vac- 
uum leads us to employ the Wilson's operator product expan- 
sion (OPE) Hi- 
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In QCD we only know how to work analytically in the per- 
turbative regime. Therefore, the perturbative part of 11(g) in 
Eq.(Q]) can be reliably calculated. However, this does not yet 
imply that all important contributions to the QCD side of the 
sum rule have been taken into account. The complete calcula- 
tion has to include the effects due to the fields of soft gluons and 
quarks populating the QCD vacuum. A practical way to calcu- 
late the vacuum-field contributions to the correlation function is 
through a generalized Wilson OPE. To apply this method to the 
correlation function ([TJ, one has to expand the product of two 
currents in a series of local operators: 

n(q) = i f dW ? *<o|rrj(*)/(0)|0> =]Tc„(e 2 )o„, (2 > 

J n 

where the set {0„} includes all local gauge invariant operators 
expressible in terms of the gluon fields and the fields of light 
quarks. Eq. (f2]) is a concise form of the Wilson OPE. The co- 
efficients C n (Q 2 ) (Q 2 = —q 2 ), by construction, include only the 
short-distance domain and can, therefore, be evaluated pertur- 
batively. Non-perturbative long-distance effects are contained 
only in the local operators. In this expasion, the operators are 
ordered according to their dimension n. The lowest-dimension 
operator with n — is the unit operator associated with the per- 
turbative contribution: C (Q 2 ) = UP er (Q 2 ), d =l. The QCD 
vacuum fields are represented in (O in the form of vacuum con- 
densates. The lowest dimension condensates are the quark con- 
densate of dimension three: Oj, = (qq), and the gluon conden- 
sate of dimension four: O4 = (g 2 G 2 ). For non exotic mesons, 
i.e. normal quark-antiquark states, the contributions of conden- 
sates with dimension higher than four are suppressed by large 
powers of K 2 q CD IQ 2 , where 1 /Aqcd is the typical long-distance 
scale. Therefore, even at intermediate values of Q 2 (~ 1 GeV 2 ), 
the expansion in Eq. (f2]i can be safely truncated after dimension 
four condensates. However, for molecular or tetraquark states, 
the mixed-condensate of dimension five: O5 = (qgcr.Gq), the 
four-quark condensate of dimension six: 0(, = (qqqq) and even 
the quark condensate times the mixed-condensate of dimension 
eight: 0% = (qqqgcr.Gq), can play an important role. The three - 
gluon condensate of dimension-six: = (g 3 G i ) can be safely 
neglected, since it is suppressed by the loop factor 1/I6n 2 . 

In the case of the (dimension six) four-quark condensate 
and the (dimension eight) quark condensate times the mixed- 
condensate, in general factorization assumption is assumed and 
their vacuum saturation values are given by: 

(qqqq) = (qqf, (qqqgcr.Gq) = (qq)(qgcr.Gq). (3) 

Their precise evaluation requires more involved analysis in- 
cluding a non-trivial choice of factorization scheme II74I1 . In 
order to account for deviations of the factorization hypothesis, 
we will use the parametrization: 

(qqqq) = P(qqf, (qqqgcr.Gq) = p(qq)(qgcr.Gq), (4) 

where p — 1 gives the vacuum saturation values and p = 2.1 
indicates the violation of the factorization assumption 117 li |75 

m. 



2.3. The phenomenological side 

The generic correlation function in Eq. (Q]i has a dispersion 
representation 



U(q 2 ) 



-I 



ds ■ 



P(s) 



r 



s + le 



(5) 



through its discontinuity, p(s), on the physical cut. The dots in 
Eq. (0) represent subtraction terms. 

The discontinuity can be written as the imaginary part of the 
correlation function: 



p(s) = -Im[Tl(s)] 
n 



(6) 



The evaluation of the spectral density (p(s)) is simpler than the 
evaluation of the correlation function itself, and the knownledge 
of p(s) allows one to recover the whole function Tl(q 2 ) through 
the integral in Eq. ([5). 

The calculation of the phenomenological side proceeds by 
inserting intermediate states for the hadron, H, of interest. 
The current j (p) is an operator that annihilates (creates) all 
hadronic states that have the same quantum numbers as j. Con- 
sequently, Tl(q) contains information about all these hadronic 
states, including the low mass hadron of interest. In order for 
the QCD sum rule technique to be useful, one must parametrize 
p(s) with a small number of parameters. The lowest resonance 
is often fairly narrow, whereas higher-mass states are broader. 
Therefore, one can parameterize the spectral density as a single 
sharp pole representing the lowest resonance of mass m, plus a 
smooth continuum representing higher mass states: 



PC*) 



A 2 5(s ■ 



m ) + Pcontis) , 



(7) 



where A gives the coupling of the current with the low mass 
hadron, H: 

(0\j\H) = A. (8) 

For simplicity, one often assumes that the continuum contri- 
bution to the spectral density, p CO nt(s) m Eq. 10, vanishes bel- 
low a certain continuum threshold so. Above this threshold, it 
is assumed to be given by the result obtained with the OPE. 
Therefore, one uses the ansatz 



Pcontis) = p (s)®(s - s Q ) . 



(9) 



2.4. Choice of currents 

Mesonic currents for open charm mesons are given in Ta- 
ble m 



Table 2: Currents for the D mesons 



state 


symbol 


current 




scalar meson 


Do 


qc 


+ 


pseudoscalar meson 


D 


iqy 5 c 


0" 


vector meson 


D* 


qjti c 


r 


axial-vector meson 


Di 




i + 



From these currents we can construct molecular currents 
which can be eigenstates of charge conjugation C and G-parity. 
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Let us consider, as an example, a current with J pc — 1 ++ for the 
molecular D°D*° system. It can be written as a combination of 
two currents 1771 17811 : 



and 



fax) = [u(x)y 5 c(x)][c(x)y, I u(x)], 



fax) = [u(x)y, 1 c(x)][c(x)y 5 u(x)]. 



(10) 



(11) 



Since the charge conjugation transformation is defined as: 
(qlc = -q T C~ 1 = q T C and (q) c = Cq T , we get 



Uu)c = -(cy 5 «)(M7/i c ) 



(fac = -{cy^iuysc) = -fa 



Therefore, the current 

jp.to 



(12) 
(13) 

(14) 



has positive C. However, this current is not a G-parity eigen- 
state. The G-parity transformation is an isospin rotation of the 
charge conjugated current: 



0> 



(fac 



-(cy 5 d)(dy^c), 
-{cjnd){dy 5 c). 



(15) 
(16) 



In the case of a charged molecular D\D* current with J p - 
0", it can also be written as a combination of two currents: 



7*1 = {cy^ysu^d-fc), 
72 = (c/M)(Jy jU y 5 c). 



(17) 
(18) 



The charge conjugation transformation in these currents leads 
to 

0*i)c = -{uy l iy5c){cy li d), 



(ji)c = -(u-fc^cy^ysd), 
and the isospin rotation gives 

0'i)g = (dy M ysc)(cy"u) = j 2 , 

{ji)G = (dY'cXcy^ysu) = ji. 
Therefore, the current 



-p O'i + h) , 
V2 



(19) 
(20) 

(21) 
(22) 

(23) 



has positive G-parity. 

In the case of tetraquark [cg][«/] currents, they can be 
constructed in terms of color anti-symmetric diquark states: 
SabAqlCVcb], where a, b, c are color indices of the SU(3) 
color group, C is the charge conjugation matrix, and T stands 
for Dirac matrices. The quantum numbers of the diquark states 
are given in Table [3] 

From these diquarks we can construct tetraquark currents 
which can be eigenstates of charge conjugation C and G-parity. 



Table 3: Currents for charmed diquark states. 



state 


current 


J p 


scalar diquark 


q'aCysCb 


+ 


pseudoscalar diquark 


q T a Cc b 


0" 


vector diquark 


qlCysy^ct 


r 


axial-vector diquark 


qlCy^Cb 


i + 



In the case of a J PC — 1 ++ current it can be written as a combi- 
nation of scalar and axial-vector diquarks: 



and 



jl = eabcededqlCysCbXqdy^Ccl), 



eabcedec(qlCyf,c b )(q d y5Ccl). 



(24) 



(25) 



It is interesting to notice that the structure of the currents in 
Eqs. d24b and d25T l relates the spin of the charm quark with the 
spin of the light quark. This is very different from the spin struc- 
ture of the heavy quark effective theory lT79ll . Heavy quark ef- 
fective theory, in leading order in 1/M, possesses a heavy-quark 
spin symmetry. Therefore, hadrons can be classified accord- 
ing to the angular momentum and parity of the light fields only. 
This gives 0~ and 1~ mesons with identical properties. 
Using the charge conjugation transformations one gets: 

Ul)c = eabcededqaysCclXqjCy^Ce) = fa (26) 

(fac = e a bcedec{q a yiiCc T b )(q T d Cy5C e ) = fa (27) 
Therefore, the current 



into 



to 



■fax))-. 



(28) 



has positive C. The i was used in Eq. d28l ) to insure that jt = j^. 
As in the case of the D°D*° molecular current, the current in 
Eq. ( 1281 is not a G-parity eigenstate. However, other combina- 
tions of tetraquark currents can be constructed to be G-parity 
eigenstates. 

In general, there is no one to one correspondence between the 
current and the state, since the current in Eq. d28l can be rewrit- 
ten in terms of a sum over molecular type currents through the 
Fierz transformation. In the appendix, we provide the general 
expressions for the Fierz transformation of the tetraquark cur- 
rents into molecular type of currents as given in Eq. (TT~4b - How- 
ever, as shown in the appendix, in the Fierz transformation of a 
tetraquark current, each molecular component contributes with 
suppression factors that originate from picking up the correct 
Dirac and color indices. This means that if the physical state is 
a molecular state, it would be best to choose a molecular type 
of current so that it has a large overlap with the physical state. 
Similarly for a tetraquark state it would be best to choose a 
tetraquark current. If the current is found to have a large over- 
lap with the physical state, the range of Borel parameters where 
the pole dominates over the continuum would be larger, and 
the OPE for the mass would have a better convergence. These 
conditions will lead to a better sum rule; this means that the 
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Borel curve has an extremum or is flat, and the calculated mass 
is close to the physical value. Therefore, if the sum rule gives 
a mass and width consistent with the physical values, we can 
infer that the physical state has a structure well represented by 
the chosen current. In this way, we can indirectly discriminate 
between the tetraquark and the molecular structures of the re- 
cently observed states. However, it is very important to notice 
that since the molecular currents, as the one in Eq. ( TPfl i. are lo- 
cal, they do not represent extended objects, with two mesons 
separated in space, but rather a very compact object with two 
singlet quark-antiquark pairs. 

When the current is fixed, we proceed by inserting it into 
Eq. (Q]i. Contracting all the quark anti-quark pairs, we can 
rewrite the correlation function in terms of the quark propa- 
gators, and then we can perform the OPE expansion of these 
propagators. For the light quarks, keeping terms which are lin- 
ear in the light quark mass m q , this expansion reads: 

s ab (x) = <orfe a (*)9*(0)]>o = T -F^- jVi 



32n 2 

8ab ._ > i<*ab ,_ 
(qq) H m a (qq)x ; 

ix 2 6, 



(■^(^ V + o*»fy ~ m q c^ v lnC-x 2 )) 
x 2 6 ab 



(qgo-.Gq) 



ab 



2 7 x3 2 



m q (qgcr.Gq)jc, 



(29) 



where we have used the fixed-point gauge. For heavy quarks, 
it is more convenient to work in the momentum space. In this 
case the expansion is given by: 

r , . t + m it^ h8 G^[^(P + m) + (P + m) ( rn 

Sab(p) - i— i — 2 6ab ~ 7 n — 2T2 

p A — m L 4 \p L — m L Y 



+ —m{g G )— r - 
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(p 2 — m 2 ) 



2\A 



(30) 



2.5. The mass sum rule 

Now one might attempt to match the two descriptions of the 
correlator: 

n" he "(Q 2 ) «-» n OPE (Q 2 ) . (3i) 

However, such a matching is not yet practical. The OPE side 
is only valid at sufficiently large spacelike Q 2 . On the other 
hand, the phenomenological description is significantly domi- 
nated by the lowest pole only for sufficiently small Q 2 , or better 
yet, timelike q 2 near the pole. To improve the overlap between 
the two sides of the sum rule, one applies the Borel transforma- 
tion 



2\n+l 



S M 2[n(q 2 )] = lim 

-q ,n—*oo 

-q 2 /n=M 2 

Two important examples are: 

!B M 2[q 2 "]=0, 

and 

&M 2 



d 

dq 2 



n(q 2 ) 



1 



(m 2 - q 2 )' 1 



1 



-nr/M 2 



(n - 1)! (M 2 ) 



2\n-\ 



(32) 



(33) 



(34) 



for n > 0. From these two results, ( f33b and ( f34t . one can see 
that the Borel transformation removes the subtraction terms in 
the dispersion relation, and exponentially suppresses the contri- 
bution from excited resonances and continuum states in the phe- 
nomenological side. In the OPE side the Borel transformation 
suppresses the contribution from higher dimension condensates 
by a factorial term. 

After making a Borel transform on both sides of the sum rule, 
and transferring the continuum contribution to the OPE side, the 
sum rule can be written as 



-m- /M 2 



dse slM 2 p OPE (s) 



(35) 



If both sides of the sum rule were calculated to arbitrary high 
accuracy, the matching would be independent of M 2 . In prac- 
tice, however, both sides are represented imperfectly. The hope 
is that there exists a range of M 2 , called Borel window, in which 
the two sides have a good overlap and information on the low- 
est resonance can be extracted. In general, to determine the 
allowed Borel window, one analyses the OPE convergence and 
the pole contribution: the minimum value of the Borel mass is 
fixed by considering the convergence of the OPE, and the max- 
imum value of the Borel mass is determined by imposing the 
condition that the pole contribution should be bigger than the 
continuum contribution. 

In order to extract the mass m without worrying about the 
value of the coupling A, it is possible to take the derivative of 
Eq. d64T i with respect to 1 /M 2 , and divide the result by Eq. d64l l. 
This gives: 



(J) 



J/° ds e--"l M2 p 0PE {s) 



(36) 



This quantity has the advantage to be less sensitive to the per- 
turbative radiative corrections than the individual sum rules. 
Therefore, we expect that our results obtained to leading order 
in a s will be quite accurate. 

2.6. Numerical inputs 

In the following sections we will present numerical results. 
In the quantitative aspect, QCDSR is not like a model which 
contains free parameters to be adjusted by fitting data. The in- 
puts for numerical evaluations are the following: z) the vac- 
uum matrix elements of composite operators involving quarks 
and gluons which appear in the operator product expansion. 
These numbers, known as condensates, contain all the non- 
perturbative component of the approach. They could in prin- 
ciple, be calculated in lattice QCD. In practice they are esti- 
mated phenomenologically. They are universal and, once ad- 
justed to fit, for example, the mass of a particle, they must have 
always that same value. They are the analogue for spectroscopy 
of the parton distribution functions in deep inelastic scatter- 
ing; ii) quark masses, which are extracted from many differ- 
ent phenomenological analyses and used in our calculation; Hi) 
the threshold parametr so is the energy (squared) which charac- 
terizes the beginning of the continuum. Typically the quantity 



yfso-m (where m is the mass of the ground state particle) is the 
energy needed to excite the particle to its first excited state with 
the same quantum numbers. This number is not well known, but 
should lie between 0.3 and 0.8 GeV. If larger deviations from 
this interval are needed, the calculation becomes less reliable. 

All in all, in QCDSR we do not have much freedom for 
choosing numbers. In the calculations discussed in the next 
sections we will use |p7ll |80[ l8ll |82|l : 



m s = (0.13 ±0.03) GeV, 
m c (m c ) = (1.23 + 0.05) GeV, 
m b (m h ) = (4.24 ± 0.06) GeV, 
(qq) = -(0.23 ± 0.03) 3 GeV 3 , 
(ss) = (0.8 ± 0.2)<^>, 



<8 2 G 2 ) 



0.88 GeV 4 



(qgcr.Gq) = m {qq), m = 0.8 GeV 



(37) 



3. The X(3872) meson 

Since its first observation in August 2003 by Belle Collab- 
oration 0]], the X(3872) represents a puzzle and, up to now, 
there is no consensus about its structure. The X(3872) has been 
confirmed by CDF, DO and BaBar JilEQ- Besides the dis- 
covery production mode B + X(3%12)K+ -> J / tC K + .the 
X(3872) has been observed in pronpt pp production B15L llol . 
However, searches in prompt e + e~ production III 811 and in e + e~ 
or yy formation lfl9ll have given negative results so far. The 
current world average mass is 



M x = (3871.20 ± 0.39) MeV 



(38) 



and the most precise measurement to date is Mx = (3871.61 ± 
0.16 ± 0.19) MeV, as can be seen by FigElQ. 



X(3872) Mass Measurements 
Belle — l 

3872.00 ± 0.60 ± 0.50 MeV/c 2 

BaBar (B + ) I ■ 

3871 .30 ± 0.60 + 0.10 MeV/c 2 

BaBar (B°) I 9 1 

3868.60 ± 1 .20 ± 0.20 MeV/c 2 

DO 1 1 

3871 .80 ± 3.1 ± 3.00 MeV/c 2 

CDF old H » 

3871 .30 ± 0.70 ± 0.40 MeV/c 2 
CDF new (preliminary) 
3871.61±0.16±0.19MeV/c 2 



old average 

3871 .20 ±0.39 MeV/c 2 

new average 

3871 .51 ±0.22 MeV/c 2 



m(D°)+m(D°*) 

3871 .81 ±0.36 MeV/c 2 



3866 



3867 



3868 3869 3870 

X(3872) Mass ( MeV/c 2 



3871 



3872 



3873 



Figure 2: An overview of the measured X(3872) masses from ref. I201I . 

The mass of the X(3872) is at the threshold for the production 
of the charmed meson pair m(D°D *) = 3871.81 ± 0.36 MeV 
1 2111 . and this state is extremely narrow: its width is less than 
2.3 MeV at 90% confidence level. 



3.1. Experiment versus theory 

Both Belle and Babar collaborations reported the radiative 
decay mode X(3872) — > yJ/tf/ 122L 12311 . which determines C 
■ Belle Collaboration reported the branching ratio: 



Eg -> J/i/ry) 

T(X -> J/lJ/7T + 7T-) 



0.14 ±0.05. 



(39) 



Recent studies from Belle and CDF that combine angular in- 
formation and kinematic properties of the n + n~ pair, strongly 
favor the quantum numbers J PC = 1 ++ I22ll24l 12511 . In particu- 
lar, in ref. 12511 it was shown that only the hypotheses J PC = 1 ++ 
and 2~ + are compatible with data. All other possible quantum 
numbers are ruled out by more than three standard deviations. 
The possibility 2~ + is disfavored by the observation of the decay 
into i//(2S )y 12611 and also by the observation of the decays into 
D°D°tt° by Belle and BaBar Collaborations 0H|. On the 
other hand, the possibility 1 ++ is disfavored by the observation 
of the decay into J/if/aj by BaBar Collaborations 08311 . In the 
following we will asume the quantum numbers of the X(3872) 
to be 1 ++ . 

From constituent quark models l84ll the masses of the possi- 
ble charmonium states with J PC = 1 ++ quantum numbers are: 
2 3 Pi(3925) and 3 3 Pi(4270), and lattice QCD calculations give 
2 3 P!(4010) B and 2 3 Pi(4067) JH. In all cases the pre- 
dictons for the mass of the 2 3 Pi charmonium state are much 
bigger than the observed mass. However, a more recent lattice 
QCD calculation gives 2 3 Pi(3853) H and, therefore, this in- 
terpretation can not be totally discarded ll88tl . In an y case, the 
strongest fact against the cc assignment for X(3872) 18418911 . is 
the fact that from the study of the dipion mass distribution in the 
X(3872) -> J/i//7T + 7r- decay, Belle lH and CDF J3 concluded 
that it proceeds through the X(3872) — > J/i[tp° decay. Since a 
charmonium state has isospin zero, it can not decay easily into 
a J/ifrp final state. 

The proximity between the X mass and the D*°D° thresh- 
old inspired the proposal that the X(3872) could be a molecular 
bound state with small binding energy 

H E III H El III 

9311 . As a matter of fact, Tornqvist, using a meson potential 
model 19411 . essentially predicted the X(3872) in 1994, since he 
found that there should be molecules near the D*D threshold 
in the J PC = 0~ + and 1 ++ channels. The only other molecular 
state that is predicted in the potential model updated by Swan- 
son S is a ++ D*D* molecule at 4013 MeV. 

In ref. J50ll Swanson proposed that the X(3872) could be 
mainly a D°D*° molecule with a small but important admix- 
ture of pJ/if/ and ajj/xfj components. With this picture the decay 
mode X(3872) — > J / 'if/ n + n" n° was predicted at a rate compa- 
rable to the X(3872) — > Jjtyn^n~ mode. Soon after this pre- 
diction, Belle Coll. ll22ll reported the observation of these two 
decay modes at a rate: 



S(X -> J/i//TT + n-n") 
S(X -» J/ifm+TT-) 



= 1.0 ±0.4 ±0.3. 



(40) 



This observation establishes strong isospin and G parity viola- 
tion and strongly favors a molecular assignment for X. How- 
ever, it still does not completely exclude a cc interpretation for 
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X since the isospin and G parity non-conservation in Eq. d40b 
could be of dynamical origin due to p° - to mixing 19511 or even 
due to final state interactions (FSI) containing D loops, such as 
X — > J /if/ to — > DD — > J/tf/p. Alth oug h all the ingredients (spe- 



cially the charm form factors H96L 19711 ) for the relevant effective 
field theory are available, there are no quantitative results in the 
FSI approach yet. 

The decay X — > J/i//a> was also observed by BaBar Collabo- 
ration 083I1 at a rate: 



S(X -> J/4m + 7T-n u ) 
S(X -> J/ifm+n-) 



0.8 ±0.3, 



(41) 



which is consistent with the result in Eq. ( l40b . 

It is also important to notice that, although a D°D* a molecule 
is not an isospin eingenstate, the ratio in Eq. d40b can not be 
reproduced by a pure D°D*° molecule. In ref. fl98ll it was shown 
that for a pure D°D*° molecule 



r(X(D u D* u ) -> Jlipn+n-ji") 
r(X(D°D*°) -> J/ipn+TT) 



= 0.15. 



(42) 



In refs. II27L 12811 Belle and BaBar Collaborations reported the 
observation of a near threshold enhancement in the D°D°7r Q sys- 
tem. The peak mass values for the two observations are in good 
agreement with each other: (3875.2 + 1.9) MeV for Belle and 
(3875.1 ± 1.2) MeV for BaBar, and are higher than in the mass 
of the X(3872) observed in the //i/ot+tT channel by (3.8 ± 1.1) 
MeV. Since this peak lies about 3 MeV above the D*°D° thresh- 
old, it is very awkward to treat it as a D*°D° bound state. Ac- 
cording to Braaten i99"ll . the larger mass of the X measured in 
the D°D°7T decay channel could be explained by the difference 
between the line shapes of the X into the two decays: D D n a 
and J/ifm + 7r~. In the decay of a narrow X molecular state into 
its constituents D*°D°, the w idth of D*° distorts the decay line 
shape of the X(3812) lIlOOll . Therefore, the peak observed in 
the B — > K D°D°n° decay channel could be a combination of a 
resonance below the D*°D° threshold from the B — > K X decay 
and a threshold enhancement above the D*°D° threshold. In this 
case, fitting the D°D°n° invariant mass to a Breit-Wigner does 
not give reliable values for the mass and width. However, in a 
new measurement 1I30I1 Belle has obtained a mass (3872.6 ±0.6) 
MeV in the D°D*° invariant mass spectrum, which is consistent 
with the current world average mass for X(3872). Using this 
new data and taking into account the universal features of the 
S -wave threshold resonance Braaten and Stapleton concluded 
that the X(387 2) is a extremely weakly-bound charm meson 
molecule [[101]. 

Other interesting possible interpration of the X(3872), first 
proposed by Maiani et al. 114511 . is that it could be a tetraquark 
state resulting from the binding of a diquark and an antidiquark 
1 95]. This construction is based in the idea that diquarks can 
form bound-states, which can be treated as confined particles, 
and us ed as degrees o f freedom in parallel with quarks then- 



selves 1102L I103L 110411 . Therefore, the tetraquark interpreta- 



tion differs from the molecular interpretation in the way that 
the quarks are organized in the state, as shown in Fig. 1. The 
drawback of the tetraquark picture is the proliferation of the 




U 





_ D antidiquark - diquark 

D°-D*° "molecule" 



Figure 3: Cartoon representation for the molecular and tetraquark interpreta- 
tions of X(3872). 



predicted states 045I1 and the lack of selection rule s that could 
explain why many of these states are not seen II 1 0511 . 

The authors of ref. fl45fl have considered diquark-antidiquark 
states with J pc = 1 ++ and symmetric spin distribution: 

X C] = [cq] s =i[cq]s=o + [cq]s=o[cq]s=i- (43) 
The most general states that can decay into 2n and 3n are: 



X\ = cos 6X U + sin 6X t j, X/, = cos 9X c i - sin 8X„. 
Imposing the rate in Eq.d40"li, they get 9 



(44) 



20°. It is important to 
notice that a similar mixture between D°D*° and D + D* molec- 
ular states with the same mixing angle 6 ~ 20° Ir98ll . would also 
reproduce the decay rate in Eq.(l40h. 

The authors of ref. 14511 also argue that if X\ dominates B + 
decays, then X/, dominates the B a decays and vice-versa. They 
have also predicted that the mass differe nce b etween the X par- 
ticle in B + and B° decays should be 1I45 II1O6I1 



M(X h ) - M(Xi) = (8 ± 3) MeV. (45) 

There are two reports from Belle l39ll and Babar Collab- 
orations for the observation of the B° — » K° X and B + — > K + X 
decays that we show in Figs. |4]and|5] 

Although the identification of the X(3872) from the decay 
B + — > K + X in these two figures is very clear, this is not the case 
for the B° — > K° X decay, where the evidence for the existence 
of such a state is not so clear. 

In any case, if one accepts the existence of the X in the decay 
B° — > K° X, these reports are not consistent with each other. 
While Belle measures 13911: 



S(B° -> XK°) 
S(B + -> XK + ) 



0.82 ± 0.22 ± 0.05, 



and 



M(X) B+ - M(X) B o = (0.18 ± 0.89 ± 0.26) MeV, 
BaBar measures [40]: 

S(B° -> XK°) 



S(B + -> XK+) 



= 0.41 ± 0.24 ± 0.05, 



(46) 



(47) 



(48) 




1.8 3.82 3.84 3.86 3.88 3.9 3.92 3.94 3.96 3.98 4 

m,, (GeV/c 2 ) 



— 80 




Figure 4: 
and B + - 



and 



3.82 3.84 3.86 3.88 3.9 3.92 3.94 3.96 3.98 4 

m,, (GeV/c 2 ) 

Fits to the m(J lipn + n~) distribution for the decays B° — » K° X (top) 
. K + X (botton) from ref.0. 



M(X) B , - M(X)bo = (2.7 + 1.6 + 0.4) MeV. 



(49) 



In both measurements, the mass difference between the two 
states is much smaller than the prediciton in Eq.d43b. 

If X(3872) is a loosely bound molecular state, the branching 
ratio for the decay B° — > K° X is suppresed by more than one 
order of magnitude co mpa red to the decay B + — > K + X. The 
prediction in refs. |6(1 107 1 gives: 



0.06 < 



F(B° -> XK°) 
T(B + XK + ) 



< 0.29, 



(50) 



which is, considering the errors, still consistent with the data 
from BaBar. 

Recentely BaBar has reported the observation of the decay 
X(3872) ip(2S )y 0] at a rate: 



S(X^^(2S)y) 



S(X ifry) 
while the prediction from ref. l50Tl gives 
r(X^if,(2S)y) 



= 3.4+1.4, 



(51) 



T(X if,y) 



4x 10" 



(52) 




r82 3.8X^.84 3.8? 3.86 3.87 3.88 3.89 
M(Jly ™t) (GeV/c 2 ) 



uj 20 




3.82 3.83 3.84 3.85 3.86 3.87 3.88 3.89 3.9 3.91 
M(JA|/ jot) (GeV/c 2 ) 

Figure 5: Fits to the m(J/i//7r + n~) distribution for the decays B° — > K° X (top) 
and B + -» K + X (botton) from ref.(Hl. 



While this difference could be interpreted as a strong point 
against the molecular model and a s a p oint in favor of a conven- 
tional charmonium interpretation [44], it can also be interpreted 
as an indication that there is a significant mixing of the cc com- 
ponent with the D a D*° molecule. As a matter of fact, the neces- 
sity of mixing a cc component with th e D°D*° molecule was al - 
ready pointed out in some works l^flMflMfTT^fTllfTlll . 
In particular, in ref. ill 1 Oh it was phenomenologically shown 
that, because of the very loose binding of the molecule, the pro- 
duction rates of a pure molecule X(3872) should be at least one 
order of magnitude smaller than what is seen experimentally. 

In ref. 15811 . a Monte Carlo simulation of the production 
of a bound D°D*° state with binding energy as small as 0.25 
MeV, obtained a cross section of about two orders of magni- 
tude smaller than the prompt production cross section for the 
X(3872) observed by the CDF Collaboration. The authors of 
ref. 115 811 concluded that S -wave resonant scattering is unlikely 
to allow the formation of a loosely bound D°D*° molecule, thus 
calling for alternative (tetraquark) explanation of CDF data. 
However, it was pointed out in ref. fl59ll that a consistent anal- 
ysis of DD* molecule production requires taking into account 
the effect of final state interactions of the D and D* mesons. 
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This observation changes the results of the Monte Carlo calcu- 
lations bringing the theoretical value of the cross section very 
close to the observed one. Thus, the question of interpretation 
of X(3872) production at hadronic machines is not y et set tled- 
Oth er interpr etations for the_XT3872) like cusp ll 1311 . hy- 
brids 1 1 14 1 15 1, or glueball 1 1 16jl have already been covered 



inrefs. [60,61,62, 63, 6J,|65| 



66L 16711 . Here we would like to 



focus on the QCD sum rules studies of this meson. 



1 ++ state we can con- 



3.2. QCDSR studies o/X(3872) 

Considering the X(3872) as a J p 
struct a current based on diquarks in the color triplet config- 
uration, with symmetric spin distribution: [cq\s=i[cq\s=o + 
l c q]s=olc<j]s="i, as proposed in ref. 14511 . Therefore, the corre- 
sponding lowest-dimension interpolating operator for describ- 
ing X q as a tetraquark state is given by: 



■(q,di) l ^abc^dec r , j„ N ,_ r^-T\ 

]J = — ■^[(q a Cy5c b ){qdy ti Cc e ) 

V2 

+ (qlCy M c b )(q d y 5 Ccl)] , 



(53) 



where q denotes a u or d quark. 

On the other hand, we can construct a current describing X q 
as a molecular DD* state: 



jf mol \x) = 



1 



(q a (x)y 5 c a (x)c b (x)y ft q h (xfj 



(q a (x)-y M c a (x)cb(x)y 5 q b (xij 



(54) 



In general, other four-quark operators with 1 ++ are possible. 
For example, starting from the simple charmed diquark states 
given in Table II, another tetraquark current with J PC = 1 ++ can 
be constructed by combining the pseudo scalar 0~ and vector 1~ 
diquark. Equivalently, additional current can be constructed for 
the meson type currents. The number of currents increases fur- 
ther, if one allows for additional color states; color sextet for the 
diquark and color octet for the molecular states. A n ext ensive 
study has been carried out for the + ^ ligh t mesons 111 1711 . with 
their mixing under renormalizations ll 1811 from which one can 
form renormalization group invariant physical currents. The 
choice of the current does not matter too much provided that 
we can work with quantities less affected by radiative correc- 
tions and where the OPE converges quite well. This is borne 
out in the well-known case of baryon sum rules, where a sim- 
ple choice of operator ll 19ll and a more general choice lll20ll 
have been studied. Even though apparently different, mainly in 
the region of convergence of the OPE, the two choices of in- 
terpolating currents have provided the same predictions for the 
proton mass. In some cases however, particular choices might 
be preferable over the others. 



X(3872). The correlator for these currents can be written as: 



IV(<7) = ij d 4 xe u '-\0\T[jl q \x)j { f\O)]\0} 

= -IW)^ - ^f) + U Q (q 2 )^, (55) 



where, since the axial vector current is not conserved, the two 
functions, III and IIo, appearing in Eq. ( IB3I ) are independent 
and have respectively the quantum numbers of the spin 1 and 
mesons. 

Using the current in Eq. (1531 1. as an example, the correlation 
function in Eq. $55[ can be written in terms of the quark propa- 
gators as: 

n„ v (g) = f d 4 xd 4 Pl d 4 P2 fi »-<™-«>{ 

2(27r) s j \ 



Tr [si ly ( Pl ) 75 CSf a ,(x)Cy 5 ] Tr [ S q , d (-x)y,CS c J e (-p 2 )Cy,] + 



+ Tr [Sl h ,( Pl )y v CSf a ,(x)C 7tl ]x 



Tr [sl d (-x)y 5 CSf e (-p 2 )C 75 ] 



(56) 



In the OPE side, we work at leading order in a s and con- 
sider the contributions of condensates up to dimension eight. 
To evaluate the correlation function in Eq. (l56*l l. we use the 
momentum-space expression for the charm-quark propagator 
given in Eq. (f30b . The light-quark part of the correlation func- 
tion is calculated in the coordinate-space, using the propagator 
given in Eq. ( |29] i. The resulting light-quark part is combined 
with the charm-quark part before it is dimensionally regular- 
ized at D = 4. 

The correlation function, III, in the OPE side can be written 
as a dispersion relation: 



d 



P(S) 



4m? s ~ q 



2 ' 



(57) 



where the spectral density is given by the imaginary part of the 
correlation functio n: np (s) = Im[Tl° PE (s)]. For the current in 
Eq. (1531 we obtain lfl2lll : 



3.2.1. Two-point function 

For the present c ase, t he tw o cur rents in Eqs. (1531 and < T54b 
were used, in refs. IU21I1 and II 12211 respectively, to study the 



p u ™(j) = p" er '{s) + p m «(s) +p {q " > (s) +p {u ~ > (s) 
+ p mi \s)+p {qq) \s) , 



(58) 
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with 



OLmax 1 —d 



1 C da T dp 

<*mm P,mn 

+P) [(a + P)m 2 c - a/3sf , 



m q f da I {qq) [m 2 - a(l - a)s] 2 



2 3 tt 4 J a \ 2 2 (I -a) 

I -a 

+ J ^[( a +flm 2 c -aps][-m 2 c {qq) + ^[(a+®m 2 c 

A™ 

a ^ + ZTTS Q + * +00 -<*-P) - apsf 

7?n L af5 l L J 



257T 4 



: [(or + P)m 2 c - Q'yS.vJ 



«„„„ A, 
2 



1—0' 



2 g 2 > r , f 48 r 2 n 



»,„,„ An 
,2n , o\2\ 



m z c (l -(a+P) 1 ) (1- 2a -2/3) 



2a 
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2 2 

(m c - a(l - a)s) 

a 



J d/3 [(a + P)nr c - aps] |^ + 



1 a+p 



(59) 



where the integration limits are given by a mi „ = 
(1 - Vl-4m?/i)/2, <w = (1 + yfl - 4m 2 /s)/2 and 
(j6m/n = am 2 ,)/(sa - m 2 ). We have included the contribution of 
the dimension-six four-quark condensate: 



(s) 



m 2 p(qq} 2 
Yin 2 



■ 4m? 



(60) 



and (for completeness) a part of the dimension-8 condensate 
contributions: 



n 



mix{qq) 
1 



(M 2 ) 



mlp{qgo-.Gq){qq) 



24n 2 



f 

Jo 



da 



1 + 



1 



2(1 -a) 





m 2 


exp 


a(l - a)M 2 



a(l - a)M 2 
(61) 



In Eqs. ( f60b and (|6TT > the parametrization in Eq. (|4]i was as- 
sumed. 

Parametrizing the coupling of the axial vector meson 1 ++ , X, 
to the current, in Eqs. (IB"3l and (l54b as: 



(0\jf\X) = A% 



(62) 



the phenomenological side of Eq 



can be written as 

- -^v + • • • , 



(63) 



where the Lorentz structure projects out the 1 ++ state. The dots 
denote higher axial-vector resonance contributions that will be 
parametrized, as explained in Sec. 2.3, through the introduction 
of a continuum threshold parameter so- After making a Borel 
transform of both sides, and transferring the continuum contri- 
bution to the OPE side, the sum rule for the axial vector meson 
X up to dimension-eight condensates can be written as: 



f/)|2 e -Mj/M 2 _ / 

J4m? 



ds e- slM ~ p 0PE {s) + U'^ q '(M z ) , (64) 



mix(qq). 



4.17 GeV 




2.0 2.2 

M 2 (GeV 2 ) 



Figure 6: The jjf ^ OPE conver genc e in the region 1.6 < M 1 < 2.8 GeV z for 



^ = 4. 17 GeV (taken from ref. 112111 '). 

For the current in Eq. 031 we show, in Fig. [6] the relative 
contribution of each term on the OPE expansion of the sum rule. 
One can see that for M 2 > 1.9 GeV 2 , the addition of a subse- 
quent term of the expansion brings the curve (representing the 
sum) closer to an asymptotic value (which was normalized to 
1). Furthermore the changes in this curve become smaller with 
increasing dimension. These are the requirements for a good 
OPE convergence and this fixes the lower limit of the Borel 
window to be M 2 > 2 GeV 2 . 

We obtain an upper limit for M 2 by imposing the constraint 
that the QCD continuum contribution should be smaller than 
the pole contribution. The maximum value of M 2 for which this 
constraint is satisfied depends on the value of sq. The compari- 
son between pole and continuum contributions for yfsj) =4.15 
GeV and p — 2.1 is shown in Fig. [7] 

Having the Borel window fixed, the mass is obtained by us- 
ing Eq. d36l l. In Fig.|8]we show the obtained X mass for different 
values of yfso and p. We can see by this figure that the effect 
of the violation of the factorization assumption, given by p, is 
similar to the effect of changing the continuum threshold. 

To study the effect of higher dimension condensates in the 
mass we show, in Fig. [9] the contributions of the individual 
condensates to Mx obtained from Eq. (l36l l. From this figure 
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2.4 2.6 

M 2 (GeV 2 ) 

Figure 7: The solid line shows the relative pole contribution (the pole contribu- 
tion divided by the total, pole plus continuum, contribution) and the dashed line 
shows the relative continuum contribution for -Jsq = 4.14 GeV and p = 2. 1. 




Figure 8: The X meson mass as a function of the Borel parameter for different 
values of -y'if) and p. The solid line and the solid line with dots show the results 
for yso =4.2 GeV using p = 1 and p = 2.1 respectively. The dashed lines 
are the same for yso = 4. 1 GeV. The crosses indicate uper limit of the allowed 
Borel window. 



we see that the results oscillate around the perturbative result, 
and that the results obtained up to dimension-5 are very close 
to the ones obtained up to dimension-8. For definiteness, the 
value of Mx obtained by including the dimension-5 mixed con- 
densate will be considered as the final prediction from the SR, 
and the effects of the higher condensates as the error due to the 
truncation of the OPE. 



The final result for Mx, obtained in 1112111 considering the 
allowed Borel window and the uncertaities in the parameters, is 



M x = (3.92 + 0.13) GeV 



(65) 



which is compatible with t he experimental value of the mass of 
the X(3872). In ref. 1 1 2 ill the mass difference in Eq.fiTb was 
also evaluated giving: 




2.3 2.4 

M 2 (GeV 2 ) 

Figure 9: The OPE convergence for Mx in the region 1.9 < M 2 < 2.8 GeV 2 for 
Sq 2 =4.1 GeV. We start with the perturbative contribution plus a very small m q 
contribution (long-dashed line) and each subsequent line represents the addition 
of one extra condensate of higher dimension in the expansion: +(qq) (solid 
line), +{g 2 G 2 ) + m q {qq) (dotted-line on top the solid line), +m 2 l (qq) (dashed 
line ), +{ qq) 2 (dot-dashed line), +m 2 l {qq) 2 (solid line with triangles) (taken from 
ref.fixita 



in agreement with the BaBar measurement in Eq. d49b . 

In the case of the current in Eq. d54b . the OPE conver- 
gence and the pole contribution determine a similar Bo rel w in- 
dow Il22ll . The result for the mass obtained in ref. IU22I1 is 
Mx = (3.87 ± 0.07) GeV, in an even better agreement with the 
experimental mass. However, due to the uncertainties inherent 
to the QCDSR method, we can not really say that the X(3872) 
is better described with a molecular type of current than with a 
tetraquark type of current. To see that we show, in Fig. [TOl the 
double ratio of the sum rules 



*mol/di — 



M 



A',;, 



(67) 



where Mx ml and Mx di are the QCDSR results obtained by using 
the currents in Eqs. ( l53l and d54l i. The ratio plotted in this figure 
was obtained by using yso =4.15 GeV and by considering the 
OPE up to dimension-6. 

We see from Fig. [10] that the deviations between the two 
QCDSR results in the allowed Borel windon are smaller than 
0.01%. 



In ref. 1112311 the same currents given in Eqs. (l53l and d54"i i 
were used to study the importance of including the width of the 
state, in a sum rule calculation. This can be done by replac- 
ing the delta function in Eq. © by the relativistic Breit- Winger 
function: 



6(s- 



m 2 ) 



1 



Fv^ 



(68) 



M(X h ) - M(X[) = (3.3 + 0.7) MeV, 



(66) 



n(s - m 2 ) 2 + sF 2 ' 
The mass and width are determined by looking at the stability 
of the obtained mass against varying the Borel parameter M 2 , 
as usual. 

Although the effect of the width is not large, in the case of 
the molecular current, Eq. d54l >, it was possible to fit the exper- 
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1.01 - 



0.99 - 
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0.95 




order to have such small decay width the authors had to make 
a bold guess about the order of magnitude of the coupling con- 
stant in the vertex XJ/if/V (where V stands for the p or u> vector 
meson): 

gx*v = 0.475. (69) 

The decay width for the decay X — > J/i[fV — » J/i //F where 
F = 7t + 7t-(k + k-7t°) for V = p(co) is given by J45i[l25ll 



2.4 2.6 

M z (GeV z ) 



Figure 10: The double ratio dmoi/di obtained from the QCDSR results for the X 
mass using the currents in Eqs. {53) and {54). 



imental mass, 3872 MeV, and the width T < 2.3 MeV simulta- 
neously with a continum threshold, yfso = 4.38 GeV, as can be 
seen by Fig. [Til 
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Figure 11: Results for the X(3872) meson mass with a molecular current. The 
crosses i ndica te lower and upper limit of the Borel window, respectively (taken 
from ref. fl23h ). 



3.2.2. Three-point function 

One important question, when proposing a multiquark struc- 
ture for the X(3872), is whether with a tetraquark or molecular 
structure for the X(3872), it is possible to explain a total width 
smaller than 2.3 MeV. As a matter of fact, a large partial de- 
cay width for the decay X — > //i/ot + 7t~ should be expected in 
this case. The initial state already contains all the four quarks 
needed for the decay, and there is no selection rules prohibiting 
the decay. Therefore, the deca y is a llowed as in the case of the 
light scalars cr and k studied in II 1 2411 . where widths of the order 
of 400 MeV were found. The decay width F(X — » J/i//tt7t), with 
X(3872) considered as a tetraquark state, was also studied in 
ref. JH]. The authors arrived at T(X -> J/i/mn) ~ 5 MeV. In 



dr 1 

— (X -» JWF) = — — 
as Xmni 



\M\ 2 B V ^ F 



Tyniv 



Pis) 



where 



Pis) 

with A(a, b, c) — a 2 + b 2 



n (s - m 2 ,) 2 + (m v T v ) 2 ' 
^A(m 2 x ,m 2 r s) 



(70) 



(71) 



2m x 

2ab - 2ac - 2bc. The invariant 



amplitude squared is given by: 

\M\ 2 = g 1 x ^ v fim x ,m f ,s), (72) 
where gxi/rv is the coupling constant in the vertex XJ/ij/V and 



f(m x , »ty, s) = - 



m 2 +s (m 2 x -m 2 ) 2 
4mt z H 



2s 



(m x - s) 



m 2 x -m 2 lJ + s 



2m 2 , 



2m\ 



(73) 



Therefore, the ratio in Eq. (l40l is given by: 

T(X J/tf/7i + n-7i°) g 2 XlP jnJ , L1 B L) ^ nnn I L) 



T{X Jlil/n+n-) 



where 



r 



ds 



{nm n y 



Pis) 

(s-m 2 v ) 2 + (m v r v ) 2 



(fim x , »v, s) 

)■ 



(74) 



(75) 



5, 



Using B^ nnn = 0.89, T w = 8.49 GeV, m M = 782.6 MeV, 
p ^ nn = 1, T p = 149.4 GeV and m p = 775.5 MeV we get 

2 

(76) 



T(X^> Jlipn + n-n°) 



T(X -> J/i//7r + n~) 



= 0.118 



The coupling constant at the XJ/if/V vertex can be evaluated 
directly from a QCD sum rule calculation, which is based on 
the three-point correlation function: 

IW(p, p', q) = j d 4 xd 4 y e ip ' * e^Tl^x, y), (77) 



with 



Tl, va (x,y) = {0\T[/jx)jV(y)%\0W), 



(78) 
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where p - p' + q and the interpolating fields are given by: 



4 _ - 

Jfi — c ay^i c a^ 

fv = ^(u a y v ii a - d a y v d a ),, 



1 



:(u a y v u a + d a y v d a ). 



With the current in Eq. d82j, the ratio between the coupling 
constants is now given by [98]: 



gx^to _ 1 ^ (cose + sine) 
gxif, p ( cos 0- sinfl)' 

(79) Using the previous result in Eq. ( |76] | we get 



(88) 



For it was shown in ref. [98] that if we use the currents in 
Eqs.(T53]l or (f54b . the QCDSR result for the ratio between the 
coupligs is given by: 



gXipu 

gX^p 



1.14. 



Using this result in Eq. ( f76*b we finally get 
T(X^ J/iffK + n-7T Q ) 



T(X J/lf/7T + 7T-) 



0.15. 



(80) 



(81) 



Therefore, to be able to reproduce the experimental result in 
Eq. d40b it is necessary to use a mixed current: 



■Ad) 



(82) 



with ja given by Eq.(l53l for a tetraquark current or by Eq.i 
for a molecular current. Using the above definitions in Eq.(f78l>. 
and considering the quarks u and d as degenerated, we arrive at 

IVC^y) = (cosfl + (-1)*' sine) Tl% a {x,y), (83) 

where N p =I p = \ and N u = 1/3, I w = 0. 

To evaluate the phenomenological side of the sum rule we 
insert, in Eq.(l78l, intermediate states for X, J /if/ and V. Using 
the definitions: 

(0\ji\J/if,(p')) = m^{p'), 
(0\jV\V(q)} = m v f v e Y (q), 
(X(p)\j x a \0) = A x e* a (p), (84) 

where A x = (cos + sin ff)A q , with A q defined in Eq. (l62l , we 
obtain the following relation: 

n <.phen)^ p , j _ '(cos 6 + sin 0)A q m^f 4 ,m v f v g Xl! , v 
(p 2 - m 2 x )(p' 2 - m 2 ){q 2 - m 2 v ) 

V mr, m, / 

where the dots stand for the contribution of all possible excited 
states, and the coupling constant, gxi//V> is defined by the matrix 
element, {J/>J/V\Xy. 

(J/^p')V(q)\X(p)} = gx^ Tallv p*e a {p)e;(p')e;{q), (86) 

which can be extracted from the effective Lagrangian that de- 
scribes the coupling between two vector mesons and one axial 
vector meson 14511 : 



T(X J I ijjn + jT- ti°) 
T(X -> J/ilfn+TT) 



a 0.15 



cose + sin( 



cos t) - sin ( 



(89) 



This is exactly the same relation obtained in refs. I45L Il25h . 
that imposes ~ 20° to reproduce the experi ment al result in 
Eq.d40li. A similar relation was obtained in ref. II 1 2711 where the 
decay of the X into two and three pions goes through a D D* 
loop. 

Using e = 20°, the XJ/i/jco coupling constant was e stimated 
from the sum rule with a tetraquark current to be IU25I1 : 



gX!/,a 



13.8 ±2.0, 



which is much bigger than the number in Eq. 
a much bigger partial decay width: 



(90) 

, and leads to 



T(X -> J/tfm + 7T-7T°) = (50 ± 15) MeV. 



(91) 



A similar width was obtained in ref. 19811 by using a molecular 
current like the one in Eq. d54l ). Therefore, from a QCDSR 
calculation it is not possible to explain the small width of the 
X(3872) if it is a pure four-quark state. 

Considering the fact that also the relation in Eq. d5TT > can 
not be reproduced with a pure molecular state, in ref. M98I1 the 
X(3872) was treated as a mixture between a cc current and a 
molecular current, similar to the mixing considered in ref. Il26ll 
to study the light scalar mesons: 



Jl{x) = sm(a)f«' mo '\x) + cos(a)/f\x), 
with jn'" 10 (x) given in Eq. ( 1541 and 



6 V2 



{qq)[c a {x)y fl y s c a {x)]. 



(92) 



(93) 



There is no problem in reproducing the experimental mass of 
the X(3872) with this current for a large range of the mixture 
angle a. Considering a in the region 5°<a<13° they get 



M x = (3.77 + 0.18) GeV, 



(94) 



which is in a good agreement with the experimental value. The 
value obtained for the mass grows with the value of the mix- 
ing angle a, but for a > 30° it reaches a stable value being 
completely determined by the molecular part of the current. 

Considering also a mixture of D + D* and D D* + compo- 
nents, the most general current is given by 



f{x) = cos 0J l '(x) + sin 0J d (x) 



(95) 



L = igx^e^id^yv^. 



(87) with J''(x) and J*(x) given by Eq. <|52}. 
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Using the current in Eq. (1951 1. the XJ/iffco coupling constant 
obtained in ref. [98] for 9 = 20° and a mixing angle in Eq. d92t 

a = 9° + 4° is: 

gx** = 5.4 ± 2.4 (96) 



which gives: 



r(x -> J/i//7r + n~n°) = (9.3 + 6.9) MeV. 



(97) 



The result in Eq. < f97b is in agreement with the experimen- 
tal upper limit. It is important to notice that the width and the 
mass grow with the mixing angle a. Therefore, there is only a 
small range for the values of this angle that can provide simul- 
taneously good agreement with the experimental values of the 
mass and the decay width, and this range is 5° < a < 13°. This 
means that the X(3872) is basically a cc state with a small, but 
fundamental, admixture of molecular DD* states. By molecu- 
lar states we mean an admixture between D°D*°, D°D*° and 
D + D* , D D* + states, as given by Eq. (1951 

In ref. jl28ll a similar mixture between a cc state and molec- 
ular states (including J/if/p and J/if/ca) was considered to study 
the X(3872) decays into J/if/y and if/(2S)y, using effective La- 
grangians. In this approach the authors only needed a small 
admixture of the cc component ( equivalent to a = 78° ± 2° 
in Eq. j92l ) to reproduce the ratio in Eq. ( BTT i. It is not clear, 
however, if with this small cc admixture it is possible to obtain 
the prompt production cross section for the X(3872) observed 
by the CDF Collaboration iBl . 

With the mixing angles a and 9 fixed, it is possible to evaluate 
the width of the radiative decay X(3872) — > J /iffy, to check if 
it is possible to reproduce the experimental result in Eq. ( 1391 . 
To do that one has just to consider the current in Eq. d95l l for 
the X(3872), and exchange the j v v current, in Eq. (l78l l. by the 
photon current j y v : 



e q qy v q , 



(98) 



q—u.d.c 



with e q = |e for quarks u and c, and e q = — \e for quark d (e is 
the modulus of the electron charge). 

The phenomenological side of the sum rule is given by 1112 



nLa"(p, p', q) 



(cos 9 + sin 9)A q m lJ ,f^e ll (p')e*(p) 
(p 2 - m 2 x )(p' 2 - m^) 



(99) 



where 



mp')\fv(q)\X(p)) = ieUq)M(X(p) -» y(q)J/i//(p')),(l00) 



with 111 1 in 

M(X(p) -» y(q)J/ifr(p')) = e e*^ (%(p)fyp')(?(q) x 

q<r 

x — (A guAgaKp ■ q + Bg„ x p K q a + Cg aK pxq^). (101) 



In Eq. dlOU . A, B and C are dimensionless couplings that are 
determined by the sum rule. Using this relation in Eq.(|99l, the 



phenomenological side of the sum rule becomes: 

nphen, , , ie(cos0 + sin6)A«»v/^ 

n^v,, (p, p,q) = —j— 2 — 2T7-72 r 

m z x (p 2 - m A x )(p' z - m^) 
X (e a ^c lcr A + ^p' A q IT q a B - £ av ^q,q a p' A C 

+ ^p\p'm c -^ p -t 

' ^' A(r p' A qAq a +p'„){A + B)^-f\. 



(102) 



The values obtained in ref. II129H for the couplings, using 9 
20° and varying a in the range 5° < a < 13" are: 



A = 18.65 + 0.94, 
A + B = -0.24 ±0.11, 

C = -0.843 + 0.008 . (103) 

The decay width is given in terms of these couplings through: 
r(X -» J/ifr y) = -'— (A + Bf + -^(A + C) 1 , 

X ijj 

where p* = (m^ - mi)/(2mx). Using the result for the decay 
width of the channel //i/ot + 7t~ in Eq. < f97T >: F(X — > J/tf/ nn) = 
9.3 + 6.9 MeV, we get 



r(X -> W y) 
T(X -> J lij/ n + n-) 



0.19 ±0.13. 



(104) 



which is in complete agreement with the experimental result in 
Eq. ([39]). Therefore, from a QCDSR point of view, the X(3872) 
is a mixture of a cc state (-97%) and molecular D°D*°, D°D*° 
(-2.6%) and D + D* , D D* + (-0.4%) states. 

3.3. Summary for X(3872) 

To summarize, there is an emerging consensus that the 
X(3872) is not a pure cc state neither a pure multi-quark state. 
From the ratio in Eq. d40l ), we know that X(3872) is not an 
isospin eigenstate, therefore, it can not be a pure cc state. On 
the other hand, the binding energy, the production rates and the 
observed ratio in Eq. dSTb are not compatible with a pure molec- 
ular state. Considering all the available experimental informa- 
tion, it is very probable that the X(3872) is a admixture of a 
charmonium state with other multi-quark states: molecular or 
tetraquark states. 

3.4. Predictions for X b , X s , X s h 

It is straightforward to extend the analysis done for the 
X(3872) to the case of the bottom quark. Using the same in- 
terpolating field of Eq. ( l53l with the charm quark replaced by 
the bottom one, the analysis done for X(3872) was repeated for 
Xi, in ref. II 1 2 lh . In this case there is also a good Borel window 
and the prediction for the mass of the state that couples with a 
tetraquark (bq)(hq) with J PC — 1 ++ current is: 

M x . = (10.27 + 0.23) GeV . (105) 
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The central value in Eq. (1147b is close to the mass of T(3S), 
and appreciably below the B*B threshold at about 10.6 GeV. 
For comparison, the molecular model predicts for Xt, a mass 
which is about 50 - 60 MeV below this threshold ll60ll . while a 
relativistic quark model without explicit (bb) clu sterin g predicts 
a value of about 133 MeV below this threshold ll 1 30fl . A future 
discovery of this state, e.g. at LHCb, will certainly test the 
different theoretical models of this state and clarify, at the same 
time, the nature of the X(3872). 

In the case of X s ([cs][cs]) and X s b ([bs] [bs]), one has to re- 
place the light quarks in the currents for X and Xt, by strange 
quarks. To extract the relatively small mass- splitting , it is ap- 



propriate to use the double ratio of moments 117 ll 1 1 3 lfl : 



M\ 



(106) 



which suppresses different systematic errors and the depen- 
dence on the sum rule parameters like sq and M 2 . The result 
obtained for this ratio in ref. 112111 is: 



= 0.984 + 0.007 . 



This leads to the mass splitting: 



M x , -M x ^-(6l± 30) MeV . 



(107) 



(108) 



Similar methods used in 117 ll 1 1 3 lfl have predicted successfully 
the values of Mo, I Mo and Mb J Mb, which is not quite surpris- 
ing, as in the double ratios, all irrelevant sum rules systematics 
cancel out. 

It is in teresting to notice that the X s mass prediction from 
ref. IU21I1 is slightly smaller than the X(3872) mass, which is 
quite unusual. Such a small and negative mass-splitting is rather 
striking and needs to be checked using alternative methods. The 
(almost) degenerate value of the X and of the X s masses may 
suggest that the physically observed X(3872) state can also have 
a cess component. 

A similar analysis can be done for the Xf {[bs][bs]) giving 
| IH: 

i — M X ' 

dl = = 0.988 ±0.018, (109) 



My 



and 



M X f - M: 



x,, 



-(123 + 182) MeV 



(110) 



We expect that the X/, -family will show up at LHCb in the 
near future, which will serve as a test of this predictions. 



4. The Y(J FL = 1"") family 

4.1. Experiment versus theory 

The e + <?~ annihilation through initial state radiation (ISR) is 
a powerful tool to search for 1 states at the B-factories. The 
first state in the 1 family disco verd using this process was the 
F(4260V It was first observed in 2005, by the BaBar Collabo- 
ration llOll . in the reaction: 



with mass M = (4259 ± 10) MeV and width F = (88 + 24) MeV. 
The 7(4260) was confirmed by CLEO and Belle Collaborations 

a. 

The BaBar Collaboration also reported a nn mass distribution 
that peaks near 1 GeV, as can be seen in Fig. Q~2] 1 10], and this 
information was interpreted as being consistent with the /o(980) 
decay. 
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Figure 1 2: Dipion mass distribution for 7(4260) -* J/Wn+n- from ref. GJj]. 

The F(4260) was also observed in the B F(4260)/T 
J '/ x i'7i + n~K~ decay ll3lll . and CLEO reported two additional de- 
cay channels: J/W^n and J/^K + R- @J. 



> 
o 



u 
> 




m(7t + 7t"J/\|/)(GeV/c 2 ) 



Figure 13: The jr + rr J/m invariant mass spectrum from ref. 132 



(111) 



In an updated report [32], BaBar has confirmed the observa- 
tion of the F(4260), shown in Fig. [13] However, the new jtjt 
mass distribution shows a more complex structure, as can be 
seem in Fig.fl4l 

Repeating the same kind of analysis leading to the ob- 
servation of the F(4260) state, in the channel e + e~ — > 
7/s«T'(25)^' + 7r ~, BaBar 13311 has identified another broad peak 
at a mass around 4.32 GeV, which was confirmed by Belle jl2tl . 
Belle found that the 4r'jT + n~ enhancement observed by BaBar 
was, in fact, produced by two distinct peaks, as can be seen in 
Fig. Q3] The masses and widths obtained by Belle and BaBar 
from fits to Breit-Wigner resonant shapes are summarized in 
Table H 
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Figure 14: Dipion mass distribution for F(4260) —> J/^V^n from ref. [32 



> 

Of 

CJ 

0) 
0) 

c 5 - 
c 
LU 



— i 1 1 r~ 



~i 1 1 1 



i 1 1 r~ 




4.5 5 5.5 

M(te + je>(2S)) (GeV/c 2 ) 



Figure 15: Evidence for the vector states 7(4360) and F(4660) from ref. 112 



Table 4: Masses and widths of 7(1 ) states, measured by Belle and BaBar 



Belle 


state 


Decay mode 


M ( MeV) T( MeV) 


7(4260) 


Jl\pn + n- 


4247 ± 108 + 19+10 


7(4360) 


i//(2S)7T + n- 


4361 ±9 ±9 74 ±15 ±10 


7(4660) 


i//(2S)tt + 7t- 


4664±11±5 48 ±15 ±3 


BaBar 


state 


Decay mode 


M ( MeV) T( MeV) 


7(4260) 


Jl\pn + n- 


4259 ±6+^ 105 ±18^ 


y(4360) 


tf/(2S)7T + 7T- 


4324 ±24 172 ±33 




0.4 0.6 0.8 

M(tiV) (GeV/c 2 ) 



0.4 0.6 0.8 1 

M(jiV) (GeV/c 2 ) 



As it is evident from Fig.[l5j there is no sign of the y(4260) in 
the \jj'n + if~ mass spectrum. The nn mass distribution reported 
in ll 2II for these two resonances can be seen in Fig. [16] 

All these three states share the following properties: they are 
vector states with the photon quantum numbers, they have large 
total widths, and they have only been observed in charmonium 
decay modes. 

Since the masses of these states are higher than the D ( * ) D ( -* ) 
threshold, if they were 1 charmonium states they should de- 
cay mainly to D^D^K However, the observed Y states do not 
match the peaks in e + e~ — > cross sections measured 

by Belle H and BaBar SQ. Besides, the ¥(3S), ¥(2£>) 
and ¥(4S ) cc states have been assigned to the well stablished 
¥(4040), ¥(4160), and ¥(4415) mesons respectively. The 
predictions from quark models for the ¥(3D) and ¥(55 ) char- 
monium states are 4.52 GeV and 4.76 GeV respectively. There- 
fore, the masses and widths of these three new Y state s are 
not consistent with any of the 1 cc states 11631 l64l 1 13211 . and 
their discovery represents an overpopulation of the charmonium 
states, as it can be seen in Fig.fPTl 

An interesting interpretation is that the y(4260) is a charmo- 
nium hybrid. Hybrids are hadrons in which the gluonic degree 
of freedom has been excited. The nature of this gluonic excita- 
tion is not well understood, and has been described by various 
models. The spectrum of char moniu m hybrids has been calcu- 
lated using lattice gauge theory 1 13311 . The result for the mass is 
approximately 4200 MeV, which is consistent with the flux tube 



Figure 16: Dipion mass distribution for the vector states F(4360) (a) and 
7(4660) (b) from ref. 0. 



mode l pred ictions II 1 3-411 . However, more recent l attice simula- 
tions 1 13511 and QCD string models calculations IU36I1 . predict 
that the lightest charmonium hybrid has a mass of 4400 MeV, 
which is closer to the mass of the y(4360). In any case, a predic- 
tion of the hybrid hypothesis is that the dominant open charm 
decay mode would be a meson pair with one S -wave D me son 
(D, D\ D s , D*) and one P-wave D meson (D u D sl ) Il37ll . In 
the case of the y(4260) this suggests dominance of the decay 
mode DD\ . Therefore, a large DD\ signal could be understood 
as a strong evidence in favor of the hybrid interpretation for the 
T(4260). In the case of the y(4360) and y(4660), since their 
masses are well above the DD\ threshold, their decay into DD\ 
should be very strong if they were charmonium hybrids. 

A critical information for understanding the structure of these 
states is wether the pion pair comes from a resonance state. 
From the di-pion invariant mass spectra shown in Figs. [14] and 
[161 there is some indication that only the y(4660) has a well 
defined intermediate state consistent with /o(980) II 1 3 8f1 . Due 
to this fact and the proximity of the mass of the if/' - £((980) 
system with the mass of the y(4660) state, in ref. 113911 . the 
y(4660) was considered as a /o(980) i//' bound state. If this 
interpretation of the y(4660) is correct, heavy quark spin sym- 
metry implies that there should be a rf c - /o(980) bound state 
114011 . This state would decay mainly into j]' c tttt, and the authors 
of ref. ll 140h predicted the mass of such state to be 4616^ MeV. 
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Figure 17: Masses of the 1 
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cc states from quark models and the masses of 



The observation of this new state would clearly determine the 
structure of the 7( 4660) . The 7(4660) was also suggested to be 
a baryonium state 114111 . a canonical 5 3 Si cc state II 14211 . and a 
tetraquark with a s c- sca lar- diquark and a sc-scalar-antidiquark 
in a 2P-wave state 1 14311 . 

In the case of 7(4260), in ref. ll46Tl it was considered as a 
sc-scalar-diquark sc-scalar-antidiquark in a F-wave state. Ma- 
iani et al. 14611 tried different ways to determine the orbital 
term and they arrived at M — (4330 ± 70) MeV, which is more 
consistent with 7(4360). However, from the nn mass distri- 
bution in refs. 11321 Il2ll . none of these two states, 7(4260) and 
7(4360) has a decay with an intermediate state consistent with 
/o(980) and, therefore, it is not clear thatjiiey should have an ss 
pair in their structure. Besides, in ref. II 1 4311 the authors show 
that the mass of a [ic]s = o[5c],5 = o tetraquark in a P-wave state 
would be 20 MeV higher than the 7(4260) mass. The authors 
of ref. 114311 . found that a more natural interpretation for the 
7(4260) would be a [qc]s=o[qc]s=o tetraquark in a f-wave state. 
The 7 (4260) was also interpreted as a baryonium A f - A f state 
1 144], a S -wave threshold effect lll45ll . as a resonance due to 
the in teraction between the three, Jlx j/nn and J/tfrKK, mesons 
jl46ll and as a 45 charmonium state 1 147 ] . Although the re ar e 
some arguments against the molecular interpretation Sim, 
the 7(4260) was also considered as a molecular state bound by 
meson exchange 1 149l 150L 151 1. The three 7 states were also 
interpreted as non-resonant manifestations of the Regee zeros 



4.2. QCDSR studies for the Y(J PL = 1~ ) states 

The Y(J PC = 1 ) states can be described by m olecular or 
tetraquark currents, with or without a ss pair. In refs. Il53lll54ll 
a QCD sum rule calculation was performed using these kind of 
currents. 



The lowest-dimension interpolating operator to describe 
a J PC = 1 state with the symmetric spin distribution: 
[C5] s=0 [c5] s= i + [ci] s= i[cs] s=0 is given by: 



€ abc-£dec r , j„ 7\ 

Jn = ^\{s a Cy 5 c b ){s d yyy 5 Cc e ) 

V2 

+ {s T a Cy5y,xC b ){s^y 5 Cc T e )] . 



(112) 




Figure 18: The j„ OPE conver gence in the region 2.8 < M 2 < 4.5 GeV 2 for 
^ = 5.1 GeV (taken from ref. fl53la 

From Fig.[l8]we see that the OPE convergence for M 2 > 3.2 
GeV 2 is better than the one shown in Fig. [6] since the perturba- 
tive contribution is the dominant one in the whole Borel region. 
Using the dominance of the pole contribution to fix the upper 
value of the Borel window, the result for the mass of the state 
described by the current in Eq 



1 12b obtained in ref. 115311 is: 

my — (4.65 + 0.10) GeV, (113) 

in excellent agreement with the mas s of the 7(4660) me- 
son. Therefore, the conclusion in ref. 115311 is that the meson 
7(4660) can be described with a diquark-antidiquark tetraquark 
current with a spin configuration given by scalar and vector 
diquarks. The quark content in the current in Eq. (Ill 2b is 
also consistent with the di-pion invariant mass spectra shown 
in Fig. Q2] which shows that there is some indication that the 
7(4660) has a well defined di-pion intermediate state consis- 
tent with /o(980). It is also very interesting to notice that the 
/o(98 0) meson may be itself considered as a tetraquark state 
1 15511 . This aspect should play an important role in the 7(4660) 
decay since, as shown in ref. I125I1 . it is very hard to explain a 
small decay width when the intital four-quark state decays into 
two final two-quark states. 

Replacing the strange quarks in Eq. ( 11 12b by a generic light 
quark q the mass obtained for a 1 state described with the 
symm etric spin distribution: [cq]s=o[cq]s=i + [cq]s=i[cq]s=o is 
I lH: 

m Y = (4.49 ±0.11) GeV , (114) 

which is bigger than the 7(4350) mass, but is consistent with it 
considering the uncertainty. 
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The 7 mesons can also be described by molecular-type cur- 
rents. In particular a D j0 (2317)D*(21 10) molecule with J PC = 
1 , could also decay into ifr'n + 7T~ with a dipion mass spectra 
consistent with /o(980). A current with J PC = 1 and a sym- 
metrical combination of scalar and vector mesons is given by: 



1 

in = -^[(s a 7^c a )(cbSb) + (c a yf,s a )(s b c b )] 
V2 



(115) 



The mass obtained in ref. 115311 for the current in Eq. dl 151 ) is 



in<p. =(4.42 ±0.10) GeV, 



(116) 



which is more in agreement with the 7(4350) mass than with 
the 7(4660) mass. 

To consider a molecular DqD* current with J pc = 1 , one 
has only to change the strange quarks in Eq. dl 151 ) by a gen eric 
light quark q. The mass obtained with such current is [ 153] 



i DoB , =(4.27 ±0.10) GeV, 



(117) 



in excellent agreement with the mass of the meson 7(4260). 
Again, in order to conclude if we can associate this molecular 
state with the meson 7(4260) we need a better understanding 
of the di-pion invariant mass spectra for the pions in the de- 
cay 7(4260) — > Jlilm + ji~ . From the spectra given in Fig.[14l it 
seems that the 7(4260) is consistent with a non-strange molecu- 
lar state D Q D*. Using a D mass ^ m Do = 2352±50 MeV, the 
DqD* threshold is around 4360 MeV and it is 100 MeV above 
the mass in Eq. ( 11 171 ). indicating the possibility of a bound state. 

A J PC =1 molecular current can also be constructed with 
pseudoscalar and axi al-ve ctor mesons. A molecular DD\ cur- 
rent was used in ref. II 5411 . The mass obtained with this current 
is: 

m D6] = (4.19 ±0.22) GeV . (118) 

Therefore, considering the errors, the molecular DD\ assigne- 
ment for the meson 7(4260) is also possible, in agreement with 
the findings of ref. 1 149], where a meson exchange model was 
used to study the 7(4260) meson, the DD\ threshold is around 
4285 MeV and very close to the 7(4260) mass, indicating the 
possibility of a loosely bound molecular state. 

4.3. Summary for Y(J PC — 1 ) states 

To summarize, the discovery of the 7(4260), 7(4360) and 
7(4660) appears to represent an overpopulation of the expected 
charmonium 1 states. The absence of open charm produc- 
tion is also inconsistent with a conventional cc explanation. 
Possible explanations for these states include charmonium hy- 
brid, tetraquark state and DqD* or DD\ molecular state for 
7(4260). The 7(4360) could be a charmonium hybrid, or a 
tetraquark state (with two axial [cs] diquarks in P-wave or with 
two scalar [cs] diquarks in f-wave). The 7(4660) could be a 
canonical 5 3 Si cc state, or a tetraquark state (with symmetrical 
[cs] s= i [cs] s=0 or with [cs] s= o[cs] s=0 in a 2P-wave). 

From the QCDSR results described is this Section, the 
7(4260) could be a DqD* or a DD\ molecular state and the 
7(4660) could be a tetraquark state with symmetrical spin dis- 
tribution: [ci]s=i[cS]5 = o + [ci]5 = o[ci]s = i. It is not possible 



to describe the 7(4360) neither as a D s qD* molecular state, 
nor as a tetraquark state with symmetrical spin distribution: 
[cq]s=\[cq] s =o + [cqh=o[cqh=i- 

5. The Z + (4430) meson 

All states discussed so far are electrically neutral. The real 
turning point in the discussion about the structure of the new 
observed charmonium states was the observation by Belle Col- 
laboration of a charged state decaying into produced in 
B + -» Kt//'n + ifll . 



5.1. Experiment versus theory 

The measured mass and width of this state are M — (4433 ± 
4 ± 2) MeV and F = (45+}*^°) MeV fll. The B meson decay 
rate to this state is similar to that for the decays to the X(3872) 
and 7(3930) mesons. There are no reports of a Z + signal in the 
J/ilm + decay channel. Since the minimal quark content of this 
state is ccud, this state is a prime candidate for a multiquark me- 
son. The Z + (4430) was observed in the ij/'n + channel, therefore, 
it is an isovector state with positive G-parity: I G — 1 + . 

Using the same data sample as in ref. 11311 . Belle also per- 
formed a full Dalitz plot analysis 13711 and has confirmed the ob- 
servation of the Z + (4430) signal with a 6.4cr peak significance, 
as can be seen in Fig. [19] The updated Z + (4430) parameters 



are: M = (4433 



■19 
■12-13 



) MeV and F = (109 



-86+74 
-43-56 



)MeV. 




MVV), GeV7c* 

Figure 19: Dalitz plot projection for the ifr'jr + invariant mass from ref. f37ll . The 
solid (dotted) histograms shows the fit result with (without) a single if/'n + state. 
The dashed histogram represents the background. 



Babar Collaboration [38] also searched the Z~(4430) signa- 
ture in four decay modes: B — > 4m~K, where if/ = J/tfs or iff' 
and K — or K + . No significant evidence for a signal peak 
was found in any of the processes investigated, as it can be seen 
inFig.EU] 

The Babar result gives no conclusive evidence of the 
Z+(4430) seen by Belle. 

There are many theoretical interpretation s o f the Z + (4430) 
structure lfl^fmfl^fl^fl60l fl6ll fl62L [Til Oil [l65 , 



I ^[l67lll68[ll69lll7dlT71 .1172H . Because its mass is close to 
the D*D\ threshold, Rosner 1 156ll suggested that it is an S - wave 
threshold effect, while others co nsidered it to be a strong candi - 
date for a D*D X molecular state lfl^[IM[l^[l^[l6ll[l62ll . 
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Figure 20: The //t/m (left) and ifr'n (right) invariant mass from ref. (38 



Other possible interpretatio ns are tetraquark state 1471 116311 . a 
cusp in the D*D\ channel fl 16411 . a baryonium state II 16511 . a 
radially excited cs state 116611 . or a hadro-charmonium state 
Ir66ll . The tetraquark hypothesis implies that the Z + (4430) 
will have neutral partners decaying into if/Tr/t]. Besides 
the spectroscopy, th ere are dis cuss ions about its production 
l lMl^[lMl7(ill and decavlfl71l. 



Considering the Z + (4430) as a loosely bound S -wave D*D\ 
molecular state, the allowed angular momentum and parity are 
J p = 0~, 1~, 2~, although the 2~ assignment is probably sup- 
pressed in the B + — » Z + K decay by the small phase space. 
Among the remaining possible 0~ and 1~ states, the former 
will be more stable as the later can also decay to DD\ in S- 
wave. Moreover, one expects a bigger mass for the J p — 1~ 
state as compared to a J p = 0~ state. There is also a quenched 
lattice QCD calculation that finds attrac tive i nteratiction in the 
D*D \ sys tem in the J p = 0~ channel 1 17311 . The authors of 
ref. II 17311 also find positive scattering lenght. Based on these 
findings, they conclude that although the interaction between 
the two charmed mesons is attractive in this channel, it is un- 
likely that they can form a genuine bound state right below the 
threshold. 



5.2. QCDSR calculations for Z + (4430) 

Considering Z + (4430) as a D*D\ molecule with J p — 0~, a 
possible current describing such state, considered in ref. II 15811 - 
is given by: 



J ! = [(d a y tl Ca)(c b y'Y 5 u b ) + (d„y IJ y5C a )(c b y' 2 u b )j . (119) 

This current corresponds to a symmetrical state D* + D^+D*°D^, 
and has positive G-parity, which is consistent with the observed 
decay Z + (4430) -> ifr'n + . 

The mass obtained in a QCDSR calculation using such a cur- 
rent was 1 158]: 



m D , Di = (4.40 + 0.10) GeV, (120) 
in an excellent agreement with the experimental mass. 



To check if the Z + (4430) could also be described as a 
diquark-antidiquark st ate w ith J p = ~ , the following current 
was considered in ref. 116211 : 



7o- = 



iec [(ulCy 5 c b )(d d Ccl) - (u T a Cc b )(d d y 5 Cc T e )] 



V2 



The mass obtained with this current was 116 



m z m 



(4.52 + 0.09) GeV, 



(121) 



(122) 



which is a little bigger than the experimental value [13], but 
still consistent with it, considering the uncertainties. Compar- 



ing this result with the result in Eq. ( 1120b , we see that the re- 
sult obtained using a molecular-type current is in a better agree- 
ment with the experimental value. As mentioned in Sec. II, this 
strongly suggests that the state is better explained as a molec- 
ular state than as a diquark-antidiquark state. This is explicitly 
borne out in the calculation for the coupling, A, between the 
state and the current defined in Eq. ([8]) of Sec. II. We get: 



A Z(o _ t = (3.75 ± 0.48) x 10~ 2 GeV 5 , 



A D , Dl = (5.66 + 1.26) x 10~ 2 GeV 



(123) 



(124) 



Therefore, one can conclude that the physical particle with 
J p = 0~ and quark content ccud couples with a larger strength 
with the molecular D*D\ type current than with the current in 
Eq. " 

Il62ll it was also considered as a diquark-antidiquark 
1~ and positive G parity: 



JT2B 
In ref. _ 
interpolating operator with J p 



V2 



[(ulCysc^iddy^ysCc]) 



+ (ulCysy^CbXddysCcl)] . 



(125) 



In this case th e Bor el stability obtained is worse than for the Z + 
with J p = 0~ II162I1 . and the obtained mass was 



mz = (4.84 + 0.14) GeV , 



(126) 



which is much bigger than the experimental value and bigger 
than the result obtained using the current with J p = 0~ in 
Eq. ( I122l i. From these results it is possible to conclude that, 
while it is also possible to describe the Z + (4430) as a diquark- 
antidiquark state or a molecular state with J p — 0~, the J p = 1~ 
configuration is disfavored. 

5.3. Summary for Z + (4430) 

A confirmation of the existence of the Z ± (4430) is critical 
before a complete picture can be drawn. If confirmed, the 
only open options for the Z + (4430) structure are tetraquark or 
molecule. The favored quantum numbers are J p = 0~. 

5.4. Sum rule predictions for B* B\ and D* S D\ molecules 

It is straightforward to extend the analysis done for the D*D\ 
molecule to the case of the bottom quark. Using the same inter- 
polating field of Eq. dl 191 ) with the charm quark replaced by the 
bottom one, the analysis done for Z + (4430) was repeated for Z b 
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in ref. ll58ll . The OPE convergence in this case is even better 
than the one for Z + (4430) and the predicted mass is 



m ZB , B = (10.74 ±0.12) GeV , 



(127) 



in a very good agreement with the prediction in ref . 117411 . 

In the case of the strange analogous meson Z + considered as 
a pseudoscalar D* S D\ molecule, the current is obtained by ex- 
changin g the d quark in Eq. (11 19t by the s quark. The predicted 
mass is 115811 : 



m D , D = (4.70 ± 0.06) GeV 



(128) 



which is bigger than the D* S D\ threshold ~ 4.5 GeV, indicating 
that this state is probably a very broad one and, therefore, it 
might be very dificult to be seen experimentally. 

6. The Z + (4050) and Z+(4250) states 

The Z + (4430) observation motivated studies of other B° — > 
K~n + (cc) decays. In particular, the Belle Collaboration has re- 
ported the observation of two resonance-like structures in the 
n + Xc\ mass distribution Q]. 

6.1. Experiment versus theory 

The two resonance-like structures, called Z + (4050) and 
Z2(4250), were observed in the exclusive process B° — > 
K~n + Xc\- The significance of each of the n + Xci structures ex- 
ceeds 5<x and, if they are interpreted as meson states, their min- 
imal quark content must be ccud. Since they were observed in 
the tt + Xc i channel, the only quantum numbers that are known 
about them are I G = 1". 
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Figure 21: The n + Xc\ invariant mass distribution from ref. 01- The solid 
(dashed) histograms shows the fit result with two (without any) n + Xcl reso- 
nance. The dots with error represent data. 

When fitted with two Breit-Wigner reconance amplitudes, 
the resonance parameters are: 

14!??) MeV, 



Mi = (4051 

T\ = (82!^!g) MeV, 
M 2 



(4248!™) MeVi 



-29- 35 ' 

r 2 = (iirj^-lt) MeV - 



The invariant mass distribution, where the contribution of the 
structure in the n + Xc\ channel is most clearly seen, is shown in 

Fig. ED 

Before this observation, it was suggested, in ref. lf5lll . that 
resonances decaying into^i and one or two pions are expected 
in the framework of the hadro-charmonium model. 

Due to the closeness of the Z ] h (4050) and (4250) masses to 
the £>*5*(4020) and £>iD(4285) thresholds, these states could 
also b e inte rpreted as molecular states or threshold effects. Liu 
et al. 117511 . using a meson exchange model find strong attrac- 
tion for the D*D* system with J p = + , while using a boson 
exchange model, the author of ref. Il76ll concluded that the in- 
terpretation of Z + (4050) as a D*D* molecule is not favored. 
In any case, it is very difficult to understand a bound molecu- 
lar state which mass is above the D*D* threshold. In the case 
of Z !(42 50), using a meson exchange model, it was shown in 
ref. 114911 that its interpretation as a fljfl or DqD* molecule is 
disfavored. 

6.2. QCDSR calculations 



In ref. 115411 . the QCD sum rules formalism was used to 
study the D*D* and D\D molecular states with I G J P = 1~0 + 
and respectively. The currents used in both cases are: 



jrrrr = (d a y ll c c )(c b 'fu b s ) , 



(130) 



and 



[(d a yny5Ca)(cby5Ub) + (d a y5c a )(cby,,y5Ub)\ . (131) 



The mass obtained with the current in Eq. ( 1 1 301 > is II 15411 : 

m D , D , = (4.15 ± 0.12) GeV (132) 

where the central value is around 130 MeV above the 
D*D*(4020) threshold, indicating the existence of repulsive in- 
teractions between the two D* mesons. Strong interaction ef- 
fects might lead to a repulsion that could result in a virtual state 
above the threshold. Therefore, this structure may or may not 
indicate a resonance. 



(129) 



For the current in Eq. ( 11311 ), the mass obtained is [ 154]: 

m DiD = (4.19 ± 0.22) GeV, (133) 

where the central value is around 100 MeV below the 
Z?iD(4285) threshold, and, considering the errors, consistent 
with the mass of the (4250) resonance. Therefore, in this 
case, there is an attractive interaction between the mesons D\ 
and D wh ich c an lead to a molecular state. 

In ref. jl23ll it was found that the inclusion of the width, in 
the phenomenological side of the sum rule, increases the ob- 
tained mass for molecular states. This means that the intro- 
duction of the width in the sum rule calculation, increases the 
mass of the states that couple to the D*D* and D\D molec- 
ular currents. As a result, using the current in Eq. ( 11311 ). it 
is possible to obtain a mass mq^d = 4.25 GeV with a width 
40 < F < 60 MeV, as can be seen in Fig. 1221 
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Figure 22: Results for the DiD molecule from ref. 112311 . Each panel shows a 
different choice of the continuum threshold. Upward and downward arrows in- 
dicate the region of the Borel window M^ in and M^ ax , respectively. Associated 
numbers in % denote the dimension eight condensate contribution for upward 
arrows and continuum contribution for downward ones. 



On the other hand, the mass of the D*D* molecule will be far 
from the 2^(4050) mass. Therefore, the authors of ref. 1 1 5411 
conclude that it is possible to describe the (4250) resonance 
structure as a D\D molecular state with I G J P — quantum 
numbers, and that the D*D* state is probably a virtual state that 
is not related with the 2^(4050) resonance-like structure. Con- 
sidering the fact that the D*D* threshold (4020) is so close to 
the 2^(4050) mass and that the t]' c '(3 1 Sq) mass is predicted to 
be around 4050 MeV (S, it is probable that the Z+(4050) is 
only a threshold effect 16511 . 

7. The y(3930) and y(4140) states 

7.1. Experiment versus theory 

The 7(3930) was first observed by the Belle Collaboration 
H in the decayfi -> #7(3930) -> KcoJ/ifr. It was confirmed 
by BaBar EH El in two channels B + -> K + toJ/tJ/ and B° -> 
K°coJ/if/. The measured mass from these two Collaborations 
are: (3943 ± 11) MeV from ref. H and (3919.1^ + 2) MeV 
from ref. [83], which gives an average mass of (3929 + 7) MeV. 
This state has positive C and G parities and the total width is 
(3 1 + 5) MeV II83I1 . The m mass distribution observed by 
BaBar is shown in Fig. [23] from where we also see the X(3872) 
observed in its decay into J/iJ/ai. 

Since the decay 7 — » //</w is OZI suppressed for a char- 
monium state 116711 . also for 7(39 30) it was conjectured t hat it 
could be a hybrid ccg state 113711 . a molecular state 1501 Il09[ 
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178L 11791 117711 . or a tetraquark state 145]. 
A recent acquisition to the list of peculiar states is the nar- 
row structure observed by the CDF Collaboration in the decay 
B + -> Y(4140)K + -> J/ip(pK + . The particle's signature peak 
can be seem in Fig. [24] The mass and width of this structure is 
M = (4143 + 2.9 + 1.2) MeV, T = (11.7+f-g ± 3.7) MeV H. 
Since the 7(4140) decays into two I G (J PC ) = 0"(1 ~) vector 
mesons, like the T(3930) it has positive C and G parities. The 
possible J p quantum numbers of a S -wave vector- vector sys- 
tem are + , 1 + , 2 + . However, since C = (-l) i+s , the J p = 1 + 
is forbidden for a state with L — and C = +1. Therefore, 
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Figure 23: Reconstruction of the X(3872) and 7(3930) peaks from their decays 
into J/t//a> mesons by BaBar (taken from ref. l83ll ). 



the possible quantum numbers for 7(3930) and T(4140) are 
J PC = ++ , 1~ + and 2 ++ . At these quantum numbers, 1~ + is 
not consistent with the constituent quark model and it is con- 
sidered exotic. 
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Figure 24: Reconstruction of the F(4140) peak from its decay into muons and 
K mesons by CDF. 

There are already some theoretical interpretations of this 
structure. Its interpretation as a conventional cc state is disfa- 
vored because, as pointed out by the CDF Collaboration Jit], it 
lies well above the threshold for open charm decays and, there- 
fore, a cc state with this mass would decay predominantly into 
an ope n cha rm pair with a large total width. It was also shown 
in ref. jl80ll that the 7(4140) is probably not a P-wave charmo- 
nium state: x' c j (J = 0, 1). If it were the case, the branching 
ratio of the hidden charm decay, 7(4140) — > J/tfxp, would be 
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much smaller than the experimental observation. 

In ref. Il8lll . the authors interpreted the 7(4140) as the 
molecular partner of the charmonium-like state 7(3930). They 
concluded that the 7(4140) is probably a D*D* molecular 
state with J PC = ++ or 2 ++ , while the 7(3930) is its D*D* 
molecular partner. This idea is supported by the fact that the 
mass difference between these two mesons is approximately 
the same as the mass difference between the <p and (±> mesons: 
rtiY(4uo) — wJy(3930) ~ m <p ~ m bj ~ 210 MeV. It is also inter- 
esting to notice that, if the 7(4140) and the 7(3930) mesons 
are D*D* and D*D* molecular states, the binding energies of 
these states will be approximately the same: niY(4\40) - 2m^ ~ 

'«y(3930) -Itnjy 90 MeV. However, with a meson exchange 

mechanism to bind the two charmed mesons, it seems natural 
to expect a more deeply bound system in the case that pions 
can be exchanged between the two charmed mesons, as in the 
D*D* , than when only 77 and <p mesons can be exchanged, as in 
the D* S D* system. 

In ref. jl82ll the author argue that the 7(4140) can be inter- 
preted either as a D*D* molecular state or as an exotic hybrid 
charmonium with J pc = 1 ~ + . A molecular DID* confi guration 
was also considered in refs. 



[183 



QCD sum rules, the authors of refs. 111851 118711 agree that it is 
possible to describe the 7(4140) with a molecular D*D* current 
with J PC = ++ . Using one boson exchange model the author of 
ref. jl86ll showed that the effective potencial of the D* S D* S sys- 
tem su pport s the explanation of 7(4140) as a molecular state. 



In ref. 118811 the 7(4140) was considered as a r 



V 



tetraquark state. The J PC — 1 ++ assignment reduces the cou- 
pling of the 7(4140) with the vector-vector channel and, there- 
fore, a small decay width would be possible in this case b ut no t 
for a J PC = ++ tetraquark state, as pointed out in ref. 1181 1. 
The authors of ref. 1118911 argue that the J/4"f> system has quan- 
tum numbers J PC = 1 and that the enhancement observed by 
CDF does not represent any kind of resonance. There is also 
a prediction for the radiative open charm decay of the 7(4140) 
that could test the molecular assignment of this state [ 190]. 



7.2. QCDSR calculation for 7(3930) and 7(4140) 

Considering the 7(3930) and 7(4140) as I G J PC = + ++ 
states, the possible currents that couple with a D*D* and D*D* 
molecular states are 

j q = (qay l iCo)(Cb'fqb) , and j s = (s a y^c a )(c b y Sb ) , (134) 

respectively. Thes e two currents were considered in a QCDSR 
study in ref. 1 1 8511 . Surprisingly the masses obtained were 
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Figure 25: The OPE convergence for the WD" molecular current, as a function 
of the Borel mass, for *>[sq = 4.4 GeV. We plot the relative contributions starting 
with the perturbative one and each other line represents the total OPE after 
adding of one extra condensate in the expansion. 



the terms in the OPE side of the sum rule, up to dimension- 
six. From this figure we see that only for M 2 > 3.5 GeV 2 
the contribution of the dimension-six condensate is less than 
20% of the total contribution. Therefore, the lower value of 
M 2 in the sum rule window should be M 2 iin = 3.5 GeV 2 . The 
inclusion of the dimension-eight condensate improves the OPE 
convergence, and its contribution is less than 20% of the total 
contribution for M 2 > 2.5 GeV 2 , as it can be seen in Fig. [26] 




m D , D , = (4.14 + 0.09) GeV, 



(135) 



Figure 26: The OPE convergence for the D*D* current for -y^ = 4.6 GeV. We 
plot the relative contributions starting with the perturbative one (long-dashed 
line), and each other line represents the total OPE after adding of one extra con- 
densate in the expansion: + (ss) (dashed line), + (g 2 G 2 ) (dotted line), + m^(ss) 
(dot-dashe d lin e), + (ss) 2 (line with circles), + m 2 (ss) 2 (line with squares). 
From ref. llS 



and 



m D , D . = (4.13 ±0.11) GeV. 



(136) 



There is another QC DSR c alculation for the D*D* molecular 
current with J PC = ++ 119111 . that finds a mass nio'D' = (3.91 ± 
0.11 ) Ge V, compatible with the 7(3930) state. However, in 
ref. Il9lll the authors have considered only the condensates up 
to dimension six. We show, in Fig. [25] the contribution of all 



Besides, we observe that considering only condensates up to 
dimension-six, the re is no pole dominance, in the Borel range 
considered in ref. 119111 . as it can be seen in Fig. [27] and, there- 
fore, no allowed Borel window can be found in this case. 

In the case of the D*D* S molecular current, the inclusion of 
the dimension-eight condensate almost does not change the 
OPE convergence and the value of the mass. This is the rea- 
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Figure 27: The pole and the continuum contributions for the D*D* current for 
V^o = 4.4 GeV. 



son why the results obtai ned w ith th e Z)*D* molecular current 
are the same in refs. 1 18511 and II 18711 . 

Ther efore , fro m a Q CD sum rule study, the results obtained 
in refs. Il85ll and lll87ll indicate that the 7(4140) narrow struc- 
ture observed by the CDF Collaboration in the decay B + — » 
Y{A\AQ)K + — > J/4npK + can be very well described by a scalar 
D* S D* current. The mass obtained with the D*D* scalar current, 
on the other hand, depends on the dimension of the condensates 
considered in the OPE side, showing that there is still no OPE 
convergence in the sum rule up to dimension-6 condensate. To 
test if the convergence is achieved up to dimension-8, it is im- 
portant to consider higher dimension condensates in the OPE 
side of the sum rule. 



8. The Z(3915) and Z(4350) states 

In recent communications, the Belle Collaboration has re- 
ported the observation of two narrow peaks in the two-photon 
processes yy — > a>J/if/ |Ql and yy — > ^J/if/ | TJ. ], as can be seen 
inFigs.|28]and|29] 
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Figure 28: The ojj/if/ m ass e nhancement observed by Belle Coll. in the yy — > 
coj/i// events. From ref. Il92ll . 

T hese two states were reported in the experimental reviews 
I68i ll92l[l93ill94ll . When fitted with Breit-Wigner resonance 
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Figure 29: The (pJ/tf> m ass e nhancement observed by Belle Coll. in the yy 
<pJ/if/ events. From ref. Il92ll . 



amplitudes, the resonance parameters are: 

M = (3914 + 4 + 2) MeV, 
r = (28 + 12!^) MeV, 



for X(3915), and 



M = (4350.6^ + 0.7) MeV, 



(13.3 



-5.1 
1-17.9 , 
-9.1 ± 



4.1) MeV. 



(137) 



(138) 



for X(4350). 

Since these two states decay into two vector mesons, they 
have positive C and G parities, like the 7(3930) and 7(4140). 
As a matter of fact, it is possible that the X(3915) be the same 
state as the 7(3930), observed by Belle fl and BaBar MUl 
in the decay channel B — > /if 7(3930) — > KutJ/i//, since the 
Babar mass and width for this state are very close to the result in 



Eq.C[37}: M = (3919.1^: 



:2MeVandT = (31^° + 5) MeV. 



In any case, a charmonium assignment for this state is difficult 

A 

In the case of the X(4350) its mass is much higher than the 
7(4140) mass and, as it can be seen in Fig. [29j no 7(4140) sig- 
nal is observed in the two- phot on process yy — > (f>J/i[s. This 
fact was interpreted, in ref. IU92I1 . as a point against the D * + D* ~ 
molecular picture for the 7(4140) states. As shown in ref. 118411 . 
a D* + D*~ molecular state should be seen in the two-photon pro- 
cess. 

The possible quantum numbers for a state decaying into // ifrcp 
are J PC = ++ , 1~ + and 2 ++ . At these quantum numbers, 1~ + 
is not consistent with the constituent quark model and it is con- 
sidered exotic. In ref. Il ill it was noted that the mass of the 
X(4350) is consistent with the prediction for a cscs tetra quar k 
state with J PC = 2 ++ lfl88ll and aD* + D*„ molecular state lfl95ll . 

with 



f D * n m olecular state 
However, the state considered in ref. II 1 9511 has J p — 1 



no definite charge conjugation. A molecular state with a vector 
and a scalar D s mesons with negat ive charge conjugation was 
studied by the first time in ref. II 1 5 311 . and the obtained mass was 
(4.42 + 0.10) GeV, also consistent with the Z(4350) mass, but 
with not consistent quantum numbers. A molecular state with a 
vector and a scalar D s mesons with positive charge conjugation 
can be constructed using the combination D* + D* - D'^D*^. 
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Some possible interpretati ons f or this state are: a excited P- 
wave charmonium state H" 2 119611 : a mixed charmonium-D*D* 
state lfl97ll . 



8.1. QCDSR calculation for X(4350) 

F ollow ing Belle Collaboration's suggestion ll ill , in 
ref. [198], a QCDSR calculation using a D* S D* Q current with 
J PC = 1~ + , was considered to test if the new observed reso- 
nance structure, X(4350), can be interpreted as such molecular 
state. 

A current that couples with a J PC = 1 
state is given by: 



D*D* sQ molecular 



Jn = [(s a y M C a )(c b S b ) - (Cay^SaXShCb)] 



(139) 



For this current the OPE convergence is very good and the 
dimension-8 condensate (obtained using the factorization hy- 
pothesis) is almost negligible, as can be seen by Fig. [30] 
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Figure 30: The OPE convergence for the J p 

the region 2.8 < M 2 < 4.8 GeV 2 for ^ = 5.3 GeV. We plot the relative 
contributions starting with the perturbative contribution plus de m, correction 
(long-dashed line), and each other line represents the relative contribution af- 
ter adding of one extra condensate in the expansion: + (ss) + m s (ss) (dashed 
line), + (g 2 G 2 ) (dotted line), + (sg<r.Gs) + m s (sgcr.Gs) (dot-dashed line), + 
(ss) 2 + m s (ss) 2 (line with circle s), + (ss)(sgcr.Gs) + m s (ss)(sg<T.Gs) (line with 
squares). Taken from ref. Il98ll . 

Considering condensates up to dimension-8 and keeping 
terms which are line ar in the strange quark mass m s , the mass 
obtained in ref. II 1 9811 was: 



m D , D , = (5.05 + 0.19) GeV, 



(140) 



where the uncertainty was obtained by considering the QCD 
parameters in the ranges: m c = (1.23 ± 0.05) GeV, m s - 
(0.13 ± 0.03) GeV, (qq) = -(0.23 + 0.03) 3 GeV 3 , m 2 Q = 
(0.8 ± 0.1) GeV 2 , and also by allowing a violation of the fac- 
torization hypothesis by using p = 2.1. 



The result in Eq. ( 1140b is much bigger than the mass of the 
narrow structure X(4350) observed by Belle. Therefore, the 



authors of ref. 1 19811 concluded that it is not possible to inter- 
pret the X(4350) as a D* s D* sQ molecular state with J PC = 1~ + . 
It is also interesting to notice that the mass obtained for a 
state described with a 1 , D*D* Q molecular current, given in 
Eq. dl 16b . is much smaller than the result obtained with the 
1~ + , D*D* sQ molecular current. This may be interpreted as an 
indication that it is easier to form molecular states with not ex- 
otic quantum numbers. 



9. The X(3940), Z (3930), Z(4160) and y(4008) states 

The X(3940) H and the X(4160) H were observed by the 
Belle Collaboration in the analysis of the M reco u(J/i[r) recoil 
spectrum in e + e~ — > J/tfr(cc). The X(3940) mass and width are 
(3942^+6) MeV and F = (37+ 2 *+8) MeV, and was observed to 



decay into D*D. The X(4160) was found to decay dominantly 
into D*D* with parameters given as M = (4156^ ± 15) MeV 
and F = (139*ij ± 21) MeV. The only other known charmo- 
nium sate observed in e + e~ — > J/i//X process has J — 0. This 
fact and the absence of these states in the DD decay suggests 
that these states favor J PC = 0~ + . However, these st ates a re ei- 
ther too low or too high to be the 77" or the rf" states IU99I1 . The 
X(4160) was identified as a 2 ++ state in two recent works: gen- 
erated from a relativistic four quark equations in ref. j200ll . and 
from a cou pled c hannel a pproa ch using a vector- vector interac- 
tion in ref. IU77I1 . In ref. 1120 111 two assignments were found to 
be possible: tj c {AS ) and the P-wave excited state Xco(3P). So 
far, no other serious theoretical attempts were performed to in- 
vestigate these states, and the nature of these states remains a 
puzzle. 

The Z(3930) was observed by Belle collaboration as an en- 
hancement in the yy — > DD event. The observed mass and 
width are M = 3929 ±5 + 2 MeV and F = 29 + 10 + 2 MeV. 
The measured properties are consistent wit h expec tations for 
the previously unseen x' c2 charmonium state 112021 18411 . 

The T(4008) was seen by the Belle collaboration as a by- 
product while confirming the 7(4260) in a two resonance fit to 
the e + e~ — > 7r + 7r~//i/( reaction Ja). There are theoretical specu- 
lation s that this state might be the t(r(3S ) or is a D*D* molecular 
state l20Hfl76ll . However, BaBar could not so far confirm ex- 
perimentally its existence |32 ] . 



10. Other multiquark states 

70.7. A D S D* molecular state 

If the mesons X(3872), Z + (4430), 7(4260) and Z+(4250) are 
really molecular states, then many other molecules should ex- 
ist. A systematic study of these molecular states and their ex- 
perimental observation would confirm its structure and provide 
a new testing ground for QCD within multiquark configura- 
tions. In this context, a natural extension would be to probe 
the strangeness sector. In particular, in analogy with the meson 
X(3872), a D S D* molecule with J p = 1 + could be formed in the 
B meson decay B — > nX s — > niJ/if/Kn). Since it would decay 
into J/iftK* — > J/iftKn, it could be easily reconstructed. 
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In ref. II 1 2211 the QCD sum rules approach was used to pre- 
dict the mass of the D S D* molecular state. Such prediction is of 
particular importance for new upcoming experiments which can 
investigate with much higher precision the charmonium energy 
regime, like the PANDA experiment at the antiproton-proton 
facility at FAIR, or a possible Super-B factory experiment. Es- 
pecially PANDA can do a careful scan of the various thresholds 
being present, in addition to going through the exact form of 
the resonance curve. 



The current used in ref. 112211 is very similar to the current in 
Eq. (O for the X(3872): 



V2 L 



Cbj5 



d b )\ 



The mass obtained for this state is 112211 : 



M DsD . = (3.96 + 0.10) GeV, 



(141) 



(142) 



which is around 100 MeV bigger than the mass of the X(3872) 
meson and below the D S D* threshold (M(D S D*) 3980 MeV). 




Figure 3 1 : The D S D* meson mass as a function of the sum rule parameter (M 2 ) 
for = 4.4 GeV (dashed line), t/sq = 4.5 GeV (solid line) and yfso = 4.6 
GeV (dot-dashed line). The crosses indicate the upper limit in the Borel region 
allowed by the dominance of the QCD pole contribution (taken from ref . 1 122]). 

The Borel curve for the mass of such state is quite stable, as 
can be seen in Fig. [31] and has a minimum within the relevant 
Borel window. Such a stable Borel curve strongly suggests the 
possibility of the existence of a D S D* molecular state with J p = 
1 + . 

Of course if such state exists, there will be no reason why 
its isospin partners would not exist as well. Therefore, charged 
states with the quark content ccsu and ecus, and with a mass 
given in Eq. ( 11421 ), should also be observed in the decay chan- 
nels: J/if/K* + mAJIxj/K*-. 

10.2. A [cc][ud] state 

Considering that the double-charmonium production was al- 
ready observed in the reaction 



it seems that it would be possible to form the tetraquark 
[cc][mJ]. Such state with quantum numbers / = 0, / = 1 and 
P = +1 which, following ref.|204], we call T cc , is especially in- 
teresting. Although the process in Eq. d 14-3b will create two cc 
pairs, the production of T cc will further involve combining two 
light anti-quarks with cc and might be suppressed. On the other 
hand, heavy ion collisions at LHC mi ght p rovide another op- 
portunity f or its non-t rivial production 120511 . As already noted 
previously 1 204 , 206ll . the T cc state cannot decay strongly or 
electromagnetically into two D mesons in the S wave due to 
angular momentum conservation nor in P wave due to parity 
conservation. If its mass is below the DD* threshold, this decay 
is also forbidden, and this state would be very narrow. 

Considering T cc as an axial diquark-antidiquark state, a pos- 
sible current describing it is given by: 



in = ^abc£d e c[{c T a Cy^c b )(.u d ysCd T e ) 



(144) 



This current represents well the most attractive configuration 
expected with two heavy quarks. This is so because the most 
attractive light antidiquark is expected to be the in the color 
tripl et, flavor anti-symmetric and spin channel 112071 1208 . 
120911 . This is also expected quite naturally from the color mag- 
netic interaction, which can be phenomenologically parameter- 
ized as, 



Vn 



c. 



-Aj ■ Ajo-j ■ ctj. 



(145) 



Here, m, A, cr are the mass, color and spin of the constituent 
quark j. Depending on the color state, the color factor A" A" is 

— I , |, — y for quark quark in the color state 3, 6 and for quark- 
antiquark in the singlet state respectively. The phenomenolog- 
ical value of Ch can be estimated from fits to the baryon and 
meson mass splitting, from which one finds that two constants 
Cm and Cg can fit all the mass splitting within the mesons and 
the baryons respectively. Al so, one fi nds that Cm is about a 
factor of 3/2 larger than C»I205, 2 1 Oil . Therefore, favorable 



J/ifr + X(3940), 



(143) 



multiquark configuration will inevitably involve diquark con- 
figuration with spin zero in color 3: this is a scalar diquark 
with anti symmetric flavor combination due to anti-symmetry 
of the total wave function involving color, spin and flavor. The 
tetraquark picture for the light scalar nonet involves states com- 
posed of two scalar diquarks. However, here, one notes that 
the two pseudo scalar mesons will have a smaller mass and the 
tetraquark state not stable against strong decay. Possible stable 
configurations emerge if either the diquark or the anti-diquark 
in the tetraquark is composed of heavy quarks; this is so because 
the large attraction between quark-antiquark will be suppressed 
by the heavy quark mass while the strong attraction in either 
the diquark or the anti-diquark will remain. Simple estimates 
within a constituent quark model predict stable configurations 
with spin zero, where a light diquark combines with a scalar 
(cb) diquark, and with spin one, whe re a light diquark com- 
bines with spin 1 heavy diquark ll210ll . The T cc falls into the 
latter case where heavy anti-diquark is a color triplet cc with 
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spin 1, as the pair should combine antisymmetrically. Although 
the spin 1 configuration is repulsive, its strength is much smaller 
than that of the light diquark due to the heavy charm quark 
mass. Therefore a constituent quark picture for T cc would be 
a light anti-diquark in color triplet, flavor anti-symmetric and 
spin (£dec[ud75CdJ]) combined with a heavy diquark of spin 
1 (eabc[cj t CjfiCb])- The simplest choice for the current which 
has a non zero overlap with such a T cc configuration is given 
inEq. ( 11441 ). While a similar configuration T ss is also possible 
1 21111 . we believe that the repulsion in the strange diquark with 
spin 1 will be larger and hence energetically less favorable. 

There are s ome predictions for the masses of the Tqq states. 
In ref. A212I1 the authors use a color-magnetic interaction, 
with flavor symmetry breaking corrections, to study heavy 
tetraquarks. They assume that the Belle resonance, X(3872), 
is a cqcq tetraquark, and use its mass as input to determine the 
mass of other tetraquark states. T hey g et Mj cc = 3966 MeV 
and M Tbb = 10372 MeV. In ref. j204ll . the authors use one- 
gluon exchange potentials and two different spatial configu- 
rations to study the mesons T cc and Tbb- They get Mt„ = 
3876 - 3905 MeV and M Tbb = 10519 - 10651 MeV. There 
are also calcu lations using expansion i n the harmoni c osc illa- 
tor basis 1 213 1. and variational method 1 214 1. In ref. l210ll . the 
authors use Eq. (1145l l together with a simple diquark picture 
to find stable heavy tetraquark configuration with spin zero T c b 
and spin one T cc , T c b, Tbb- 

In ref. ll215ll a QCD sum rule calculation was done with the 
current in Eq. (1144b . The authors found that the results are not 
very sensitive to the value of the charm quark mass, neither to 
the value of the condensates. The QCDSR predictions for the 
T cc mesons mass is: 



M Tcc = (4.0 + 0.2) GeV, 



(146) 



in a very good agreement with the pred ictions based on the 
one gl uon e xchange potential model 120411 . and color-magnetic 
model j2l3l . 

It is straightforward to extend the analysis done for the T cc 
to the the bottom sector. The prediction for the Tbb mass from 
ref. Jill is: 

M Tbb = (10.2 ± 0.3) GeV , (147) 



also in a very g ood agreement with the results in refs. [204], 
EH and fcrf . 



The results from ref. 121511 show that while the T cc mass is 
bigger than the D*D threshold at about 3.875 GeV, the T bb 
mass is appreciably below the B*B threshold at about 10.6 GeV. 
Therefore, these results indicate that the Tbb meson should 
be stable with respect to strong interactions and must decay 
weakly. These results also confirm the naive expectation that 
the exotic states with heavy q uarks tend to be more stable than 
the corresponding light states Il216ll . 



11. Summary 

We have presented a review, from the perspective of QCD 
sum rules, of the new charmonium states recently observed by 
BaBar and Belle Collaborations. As it was seen case by case, 



this method has contributed a great deal to the understanding of 
the structure of these new states. When a state is first observed 
and its existence still needs confirmation, a QCDSR calculation 
can be very useful. It can provide evidence in favor or against 
the existence of the state. We have computed the masses of 
several X, Y and Z states and they were supported by QCDSR. In 
some cases a tetraquark configuration was favored and in some 
other cases a molecular configuration was favored. Q CDSR 
does not always corroborate previous indications. In Ref. 121711 . 
a comprehensive QCDSR analysis of light tetraquarks led to 
the conclusion that their existence (as tetraquarks states) is very 
unlikely. However, a more positive conclusion was found in the 
case of tetraquarks with at least one heavy quark. 

The limitations in statements made with QCDSR estimates 
come from uncertainties in the method. However these state- 
ments can be made progressively more precise as we know 
more experimental information about the state in question. One 
good example is the X(3872), from which, besides the mass, 
several decay modes were measured. Combining all the avail- 
able information and using QCDSR to calculate the observed 
decay widths, we were able to say that the X(3872) is a mixed 
state, where the most important component is a cc pair, which 
is mixed with a small molecular component, out of which a 
large fraction is composed by neutral D and D* mesons with 
only a tiny fraction of charged D and D* mesons. This con- 
clusion is very specific and precise and it is more elaborated 
than the other results presented addressing only the masses of 
the new charmonia. This improvement was a consequence of 
studying simultaneously the mass and the decay width. This 
type of combined calculation will eventually be extended to all 
states. 

We close this review with some conclusions from the results 
presented in the previous sections. They are contained in Tables 
Hand|6] 

In Table [5] we present a summary of the most plausible inter- 
pretations for some of the states presented in TableQ] described 
in the previous sections. 

Table 5: Structure and quantum numbers from QCDSR studies. 



state 


structure 


J™ 


X(3872) 


mixed cc - DD* 


1 ++ 


F(4140) 


D*D* molecule 


++ 


Z 2 + (4250) 


DD\ molecule 


1" 


7(4260) 


D()D* or DD i molecule 


1— 


Z + (4430) 


D*D\ molecule 


0" 


T(4660) 


[ci][cS] tetraquark 


1— 



In Table [6] we present some predictions for bottomonium 
states, that may be observed at LHCb. 

From the last Table we see that the mass predictions for the 
states Xb and Tbb are basically the same. These states, if they 
exist, may be produced at LHCb. Since the Tbb is a charged 
state, it could be more easily identified. Moreover, since its 
mass prediction lies below the BB* threshold, it should be very 
narrow because it could not decay strongly. 
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Table 6: Prediction for bottomonium states from QCDSR studies. 



state 


structure 


J™ 


mass (GeV) 


x b 


[bq] [bq] tetraquark or 
BB* molecule 


1 ++ 


10.27 ± 0.23 


z b 


[bq] [bq] tetraquark or 
B*B\ molecule 


0" 


10.74 ±0.12 


Tbb 


[bb][ud] tetraquark 
BB* molecule 


1 ++ 


10.2 ±0.3 



Looking at Table [6] we see that the predictions made with 
tetraquark and molecular currents are the same. Probably, in 
order to distinguish one from the other, we will need to know 
the decay width of these states. On the other hand, from the 
point of view of measurement, this convergence of predictions 
makes them robust. In other words, even not knowing very pre- 
cisely their structure, we know that these states must be there. 

Table |5] represents the final result of a comprehensive effort 
and a careful analysis of several theoretical possibilities in the 
light of existing data. It is an encouraging example of what 
QCDSR can do. This Table contains a short summary of what 
we have learned about the new charmonium states in the recent 
past. Table [6] addresses the near future and shows some predic- 
tions. 

From what was said above, hadron spectroscopy is and will 
be a very lively field and QCDSR is a powerful tool to study it. 



and the sum rule will not be able to reproduce the mass well. 
The opposite will also be true. 

To show the suppression, we discuss Fierz transformation of 
tetraquark current made of diquark fields into quark-antiquark 
fields. The starting relation is given as follows. 



<i% = -Ik 



ah 



-irr'HqyrfKq) + y^(qy 5 A a q) 



+ ^J(qo- f , v A a q) 
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Sijiqq) + yiiiqy^q) 



-(y tI y 5 )ij(qy fl y 5 q) + y-Xqy 5 q) 



+ ~v J *J(q°~ f iv<i) 



(148) 



Here, a,b are the color and i,j the Dirac spinor indices re- 
spectively. The two color representations come from 3 (^) 3 = 
801. ' 

In the formulas to follows, the two color terms can be worked 
out using the following formulas, 

£abcedec{q T S x 5 bd c T e ){qA a T2c) = 2{cT\ q)(qT 2 c) (149) 
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12. Appendix: Fierz Transformation 

In general, the currents constructed from diquark antidiquark 
type of fields are related to those composed of meson type of 
fields. However, the relation are suppressed by a typical color 
and Dirac factors so that we obtain a reliable sum rule only 
if we have chosen the current well to have a maximum overlap 
with the physical meson. This is expected to be particularly true 
for multiquark configuration with special diquark or molecular 
structures. Therefore, if we obtain a sum rule that reproduces 
the physical mass well, we can infer the inner structure of the 
multiquark configuration. For example, if the state is of molec- 
ular type, it will have a maximum overlap with a current that is 
constructed with the two corresponding meson current; the sum 
rule will then be able to reproduce the physical mass well. On 
the other hand, if a current with a diquark-antidiquark type of 
current is used, the overlap to the physical meson will be small 



We start by applying the formula to the current composed of 
two scalar diquarks. After Fierz transformation, we find, 



eabcedec{q T a Cysc b ){qdysCc T e ) 



- (cA a q)(qA a c) - (cA a y,q){qA a y,c) 



-(cA a y fl y 5 q)(qA a y f ,y 5 c) 
-(cA a y 5 q)(qA a y 5 c) + ^(cA a o- flv q)(qA a o-, J yC) 



4 



(cq)(qc) - (cy^qXqy^c) 



-{cy^ysq^qy^c) 
-(cy 5 q)(qy 5 c) + ^ica-^q^qcr^c) 



(150) 



This is a typical example of the expansion. The factor of 1/6 
multiplying the color singlet quark-antiquark types of currents, 
comes from the color and Dirac factors and is responsible for 
the small overlap to the various quark-antiquark type of cur- 
rents. In Eq. ( 1150b . all meson type of currents contributes. 
In general however, depending on the current type and charge 
conjugation, only certain types of meson-meson currents con- 
tributes. Several examples for the expansion used in the text are 
given below. 
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1. ForX(3872) with J = 1 ++ , we find, 

■(q,di) t^abc^dec r , T r< \/- r*-T\ 

$ = ^[(q a Cy5C b )(q d y^Cc e ) 

V2 

+{q T a Cy lJl c b ){q d ysCc T e )] 

= ^f ( 4 )( " i)2 [x(^ crrcr5 

+ r5 cr r c^)?)(^'Tc)] 
+ r5 cr r c r/J ) ? ^rc) 



Jft 



i 1 



2(cA a y 5 q)(qA a y,c) 
-2(cA a yil q){qA a y 5 c) - liicA^^q) 
x(qA a y v y 5 c) + 2i(cA a y v y i q)(qA a <T ll yC) 



V2 6 



2{cy 5 q)(qy M c) 
-2(cy^q)(qy 5 c) - 2i{ccr llv q)(qy v y 5 c) 



+2i(cy v y 5 q)(qcr llv c) 



(151) 



numerically, the ratio A^/Ao'Dt = 0.66. On the 
other hand, assuming that the Z + (4430) has only a 
molecular component that couples dominantly to the 
molecular current, the ratio should be 1/3, as can be 
seen from Eq. (11521 i. The fact that it is larger than 
1/3 suggest that while the Z + (4430) is dominantly of 
molecular type, there still remains some coupling to 
other type of non-molecular type of current compo- 
nent, such as the diquark-antidiquarktype. 

with J p = 1~, we find the following current does not 
fit any of the observed Z states. 



Jk = 



V2 



[(ulCysc^iddy^ysCcl) 



+{u T a Cy tl ysc b ){d d j5Cc T e )] 

- 2(cA a y 5U )(dA a y,y 5 c) 
-2{cA a y M y 5U )(dA°y 5 c) 

-2i(cA a cr flv u)(dA a y v c) - 2i(cA a y v u)(dA a <T llv c) 

- 2(cy 5 u)(dy^y 5 c) - 2(cy ll y 5 u)(dy 5 c) 



V2 6 



-2/(ccr /JV M)(c/y v c) - 2i(cy v u){d(T IIV c) 



(153) 



As can be seen from the above equation, the molecu- 



lar component as given in f^' mol) 



(cy^iqy^c) 



(cy^q)(qysc) 



comprises only a small part of the wave 



function given in Eq. ( I15U . The overlap to the molecu- 
lar type is suppressed by | so that the overall contribution 
to the correlation function is suppressed by ^. 
2. For Z + (4430) 

• with J p = ~, we find the following current works 
well. 



■/() 



(di) _ t^abc^dec 



[(«„ Cy 5 c h )(d d Cc e ) 



V2 

-(u T a Cc b ){d d ysCc T e )] 



2(cA a y fl y 5 u)(dA a y f ic) 
y M u)(dA a 7 u 75 c) 
2(cy M r5")(^7/jc) 



V2 6 



+2(cy l _ 1 u)(dy IJ y 5 c) 



(152) 



As can be seen from the above equation, the molecu- 
lar component as given in Eq. (II 19l l. contributes with 
a suppression factor of 1/3. It is useful to compare 
this factor to the ratio of the overlaps found in QCD 
sum rule analysis and given in Eqs. (1 1 23b and ( 1124b ; 
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